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ABSTRACT 


Orthogonal frequency division multiplexing (OFDM) is being successfully used 
in numerous applications. It was chosen for IEEE 802.11a wireless local area network 
(WLAN) standard, and it is being considered for the fourth-generation mobile communi¬ 
cation systems. Along with its many attractive features, OFDM has some principal draw¬ 
backs. Sensitivity to frequency errors is the most dominant of these drawbacks. In this 
thesis, the frequency offset and phase noise effects on OFDM based communication sys¬ 
tems are investigated under a variety of channel conditions covering both indoor and out¬ 
door environments. The simulation performance results of the OFDM system for these 
channels are presented. 


v 



THIS PAGE INTENTIONALLY LEFT BLANK 


vi 



TABLE OF CONTENTS 


I. INTRODUCTION.1 

A. OBJECTIVE.1 

B. RELATED RESEARCH.2 

C. ORGANIZATION OF THE THESIS.2 

n. OFDM-BASED IEEE 802.11A STANDARD.3 

A. INTRODUCTION TO OFDM.3 

B. FUNDAMENTALS OF OFDM.6 

1. Multi-Path (Delay-Spread or Time Dispersion).6 

2. Orthogonality.9 

3. OFDM Transmitter..10 

a. Channel Coding . 11 

b. Interleaving . 11 

c. Symbol Mapping . 12 

d. Differential Quadrature Phase Shift Keying (DQPSK) 

Modulation . 13 

e. IFFT . 13 

f Guard Interval Addition . 14 

g. Number of Sub-carriers . 16 

h. Symbol Pulse Shaping (Windowing) . 17 

i. RF Modulation . 17 

4. OFDM Receiver..17 

a. Removing Guard Interval and FFT Processing . 18 

b. Decoding and Deinterleaving . 18 

c. Received message and Viterbi Decoding . 18 

C. SUMMARY.19 

III. FREQUENCY ERRORS IN OFDM AND THEIR EFFECTS.21 

A. FREQUENCY ERRORS.21 

B. CREATING A SIMPLE REPRESENTATION OF AN OFDM 

SIGNAL TO SIMPLIFY ANALYSIS.23 

C. THE EFFECT OF DOPPLER FREQUENCY SHIFT, PHASE 
NOISE AND OSCILLATOR FREQUENCY OFFSET ON THE 

OFDM SYSTEMS.23 

1. Effect of Doppler Shift on the Carrier Frequencies, Sub¬ 
carriers, Envelope and Symbol Timing.24 

2. Phase Noise, Oscillator Frequency Offset and Their Effects.25 

a. Phase Noise .27 

b. Frequency Offset . 28 

D. SUMMARY.30 

IV. OFDM SYSTEM SIMULATIONS IN MATLAB.31 

A. SELECTING SYSTEM PARAMETERS.31 

vii 






































B. SIMULATION METHODOLO GY.32 

1. Code Check (Channel Model 0).34 

2. OFDM Performance in an AWGN Channel (Channel Model 1) ..35 

3. OFDM Performance in Mobile Channel (Channel Model 2).37 

a. Mobile Channels . 37 

b. Mobile Channel Simulations: . 38 

4. Performance with Frequency offset.43 

a. Noise Free Channel . 43 

b. AWGN Channel . 46 

c. Mobile Channel . 47 

5. Performance with Phase Noise.48 

a. Noise-Free Channel . 48 

b. AWGN Channel . 51 

c. Mobile Channel . 52 

C. SUMMARY.53 

V. CONCLUSION.55 

A. SUMMARY OF THE WORK DONE.55 

B. SIGNIFICANT RESULTS AND CONCLUSIONS.55 

C. SUGGESTIONS FOR FUTURE STUDIES.56 

APPENDIX A. MATLAB CODE EXPLANATION.59 

A. TRANSMITTER.59 

1. CDRCDLFT.59 

2. TDA.61 

B. CHANNEL.61 

1. Phase Noise Simulations.61 

2. Channel Model Simulations.62 

a. Channel Model 0 . 62 

b. Channel Model 1 . 62 

c. Mobile Channel Models . 63 

3. Frequency Offset Simulations.65 

C. RECEIVER.65 

APPENDIX B. MATLAB CODE.67 

LIST OF REFERENCES.127 

INITIAL DISTRIBUTION LIST.131 


viii 




































LIST OF FIGURES 


Figure 1. A Multi-Path Scenario.6 

Figure 2. Delayed Signals (After Ref. 13.).6 

Figure 3. Representation of a Symbol in A Frequency Selective Channel: (a) Time 

domain (b) Frequency domain.7 

Figure 4. Illustration of ISI (After Ref. 14.).7 

Figure 5. (a) Time Domain Representation, (b) Frequency Domain Representation 

of a Symbol in Flat Fading Channel.8 

Figure 6. OFDM Splits a Data Stream into N Parallel Data Streams (After Ref. 15.).8 

Figure 7. Representation of an OFDM Sub-Carrier: (a) In Time Domain (B) 

Frequency Domain.9 

Figure 8. Sub-Carrier Orthogonality (After Ref. 16.).10 

Figure 9. Schematic Diagram of an OFDM Transmitter.10 

Figure 10. Convolutional Encoder with Constraint Length= 7. Each box represents a 

shift register (After Ref. 11.).11 

Figure 11. 16-QAM and QPSK Constellations (From Ref. 11.).12 

Figure 12. Effects of Guard Interval (After Ref. 16.).16 

Figure 13. Schematic Diagram of an OFDM Receiver.17 

Figure 14. OFDM Receiver Front End (After Ref. 3.).22 

Figure 15. Effects of Frequency Offset: Reduced Amplitude and ICI (After Ref. 24.) ....25 

Figure 16. SNR Degradation Versus - 3-dB Bandwidth of the Phase Noise Spectrum 

for QPSK.28 

Figure 17. SNR Degradation Versus the Normalized Frequency Offset for QPSK.29 

Figure 18. Matlab Simulations Scheme.31 

Figure 19. (a) Transmitted and (b) Received QPSK Constellations After Differential 

Decoding, from a Channel Model 0 Simulation.34 

Figure 20. Effect of AWGN on Received Signal Constellation for a = 0.02 (Before 

Differential Decoding).36 

Figure 21. BER versus E b /N 0 Plot for a QPSK Simulation in AWGN Channel.37 

Figure 22. Modeling of Mobile Channels 1 and 2 Using FIR Filters and Delay 

Elements D.38 

Figure 23. Effect of Mobile Channels on Received Signal Constellations (Before 
Differential Decoding), (a) Mobile Channel 1 (b) Mobile Channel 2 (c) 

Mobile Channel 3 (d) Mobile Channel 4.40 

Figure 24. Received Signal Constellation Affected by Mobile Channels (After 

Differential Decoding), (a) Mobile Channel 1 (b) Mobile Channel 2 (c) 

Mobile Channel 3 (d) Mobile Channel 4.41 

Figure 25. Effect of Mobile Channels on the Magnitude Plot of Received Symbols. 

(a) Mobile Channel 1 (b) Mobile Channel 2 (c) Mobile Channel 3 (d) 

Mobile Channel 4.42 

Figure 26. BER versus E b /N 0 Plots for QPSK Simulation in Mobile Channels.43 


IX 




























Figure 27. Received Signal Constellation in Channel Model 0 with Frequency Offset 

of 0.1% of Frequency Spacing (Before Differential Decoding).45 

Figure 28. Received Signal Constellation in Channel Model 0 with Frequency Offset 

of 0.1% of Frequency Spacing (After Differential Decoding).45 

Figure 29. BER versus E b /N a Plots for QPSK in AWGN Channel with Relative 

Frequency Offset Values from 0% to 0.4%.46 

Figure 30. BER versus E b /N a Plots for QPSK in Mobile Channel 1 with Relative 

Frequency Offset Values from 0% to 0.4%.47 

Figure 31. Phase Noise Generated as a Wiener Process.49 

Figure 32. Received Signal Constellation in Channel Model 0 with p / R = 0.0001. 

(Before Differential Decoding).50 

Figure 33. Received Signal Constellation in Channel Model 0 with p / R = 0.0001. 

(After Differential Decoding).50 

Figure 34. BER versus E b /N a Plots for QPSK in AWGN Channel with Norma li zed 

Bandwidth P/R s Values from 0.0% to 0.28%.51 

Figure 35. BER versus E b /N a Plots for QPSK in Mobile Channel 1 with Norma li zed 

Bandwidth P/R s Values from 0.0% to 0.28%.53 

Figure 36. Hierarchical m-file Block Diagram for the Transmitter.61 

Figure 37. Phase Noise Creating Block in the m-file Hierarchy.62 

Figure 38. Channel Model 0 m-file Hierarchy.62 

Figure 39. Mobile Channels m-file Hierarchy.64 

Figure 40. Summary of Channel Model m-file Architecture.64 

Figure 41. Frequency Offset Block in the m-file Hierarchy.65 

Figure 42. Hierarchical m-file Block Diagram for the Receiver.66 

Figure 43. General m-file Block Diagram.66 


x 




















LIST OF TABLES 


Table 1. OFDM Timing Related Parameters (From Ref. 11.).4 

Table 2. OFDM Rate Dependent Parameters (From Ref. 11.).4 

Table 3. Major Parameters of OFDM Physical Layer (From Ref. 11.).5 

Table 4. Simulation Steps.33 

Table 5. Code Check Simulation Parameters.34 

Table 6. AWGN Channel Simulation Parameters.36 

Table 7. FIR Filter Coefficients for Mobile Channels 1 and 2.38 

Table 8. Mobile Channel, Simulation 1 Parameters.39 

Table 9. Frequency Offset in Noise Free Channel Simulation.44 

Table 10. Frequency Offset in AWGN Channel Simulation.46 

Table 11. Phase Noise in Noise Free Channel Simulation.48 

Table 12. Phase Noise in AWGN Channel Simulation.51 


xi 















THIS PAGE INTENTIONALLY LEFT BLANK 



LIST OF ACRONYMS 


ADSL 

AWGN 

BER 

CDMA 

DAB 

dB 

DFT 

DSSS 

DVB-T 

FCC 

FDMA 

FEC 

FHSS 

GSTN 

HDSF 

HDTV-T 

ICI 

IDFT 

IF 

IFFT 

ISI 

MAC 

MC-CDMA 

OFDM 

PSK 

QAM 

QPSK 

RF 

TDM A 
UMTS 
VHDSF 
WEAN 


Asymmetric Digital Subscriber L ine 
Additive White Gaussian Noise 
Bit Error Rate 

Code Division Multiple Access 
Digital Audio Broadcasting 
Decibel 

Discrete Fourier Transform 
Direct Sequence Spread Spectrum 
Digital Video Broadcasting-Terrestrial 
Federal Communications Commission 
Frequency Division Multiple Access 
Forward Error Correction 
Frequency Hopping Spread Spectrum 
General Switched Telephone Network 
High-Rate Digital Subscriber F ine 
High Definition TV-Terrestrial 
Inter Carrier Interference 
Inverse Discrete Fourier Transform 
Intermediate Frequency 
Inverse Fast Fourier Transform 
Inter Symbol Interference 
Medium Access Control 
Multi Carrier CDMA 

Orthogonal Frequency Division Multiplexing 

Phase Shift Keying 

Quadrature Amplitude Modulation 

Quadrature PSK 

Radio Frequency 

Time Division Multiple Access 

Universal Mobile Telecommunications System 

Very High-Rate Digital Subscriber F ine 

Wireless Local Area Network 


xiii 



THIS PAGE INTENTIONALLY LEFT BLANK 


xiv 



ACKNOWLEDGMENTS 


First of all, I want to thank my girlfriend Guzelyasemin GULAY for being in my 
life and for her support whenever I needed. I also thank my thesis advisors Prof. Murali 
Tummala and Prof. Roberto Cristi for their patience in answering my questions and for 
the enjoyable and encouraging atmosphere they created during my whole thesis process. 
Finally, I want to thank my family beginning with my father Hasan ERDOGAN, who is 
my inspiration in life, my mother Yadigar ERDOGAN, the strongest woman in my world 
and my sister llk nur ERDOGAN along with her son Mert for being a part of my life. 

I dedicate this thesis to my family. I feel really lucky to have these wonderful 
people in my life. 


xv 



THIS PAGE INTENTIONALLY LEFT BLANK 


xvi 



EXECUTIVE SUMMARY 


Due to the fact that the need to send more information increases each day, the 
military is also in need of solutions to make more information transmission to the battle¬ 
field possible. The military is looking wireless ADHOC networks as a solution to this 
problem and orthogonal frequency division multiplexing (OFDM) is the technique used 
in ADHOC networks. 

OFDM is being successfully used in numerous applications. It was chosen for 
IEE E 802.1 la wireless local area network (WLAN) standard, and it is being considered 
for the fourth-generation mobile communication systems. Along with its many attractive 
features, OFDM has some principal drawbacks. Sensitivity to frequency errors is the 
most dominant of these drawbacks. 

The objective of this thesis was to investigate the frequency offset and phase 
noise effects on OFDM based communication systems to determine the ways by which 
we can use this technique in battlefield applications more efficiently. A comprehensive 
background of OFDM based communication systems is presented, and the sub-carrier 
orthogonality requirement, inverse fast Fourier transform and fast Fourier transform 
processes are discussed. Theoretical analysis of frequency errors and their effects on the 
overall system performance are examined based on the work reported in the literature. 

Matlab simulations of an OFDM communication system were performed for the 
following channels: ideal, additive white Gaussian noise (AWGN), and mobile. The mo¬ 
bile channel selected includes multipath and Doppler effects in AWGN. The effects of 
frequency offset and phase noise were studied for all channel models by observing the 
performance deviation from the AWGN case as a function of increasing frequency offset 
and phase noise. The bit error rate (BER) performance plots of the OFDM system af¬ 
fected by frequency errors under AWGN and mobile channel conditions were obtained 
through simulations. 
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It was observed that the OFDM system’s performance degraded as a power of the 
relative frequency offset. A relative frequency offset of 0.1% required a bit energy to 
noise power spectral density ratio ( E b /N o ) value of 1 dB more than the AWGN case 

without offset, and a frequency offset of 0.2% required 2.8 dB more E b /N a for the same 
BER. 

In the phase noise case, it was observed that the OFDM performance is degraded 
linearly as we increased the normal is ed bandwidth, as was suggested in the literature. A 
norma li zed bandwidth of 0.07%, representing the amount of phase noise, required an 
E b /N 0 value of 1 dB more than the AWGN case without phase error, and each addi¬ 
tional 0.07% increment in the norma liz ed bandwidth resulted in 1 dB more E h /N 0 re¬ 
quirement. 

The effects of mobility were also studied. The mobility is simulated by including 
the Doppler frequency shift and multipath effects. The simulations suggested that the sig¬ 
nal to noise ratio degradation due to the Doppler shift, most of the time, is negligible. The 
degrading effects of multipath are observed by comparing the system performance under 
mobile channels with that under AWGN conditions. 



I. INTRODUCTION 


Orthogonal frequency division multiplexing (OFDM) is being successfully used 
in numerous applications, such as European digital audio broadcasting and digital video 
broadcasting systems [1, 2]. In 1999, the IEEE 802.11a working group chose OFDM for 
their 5-GHz band wireless local area network (WLAN) standard, which supports a vari¬ 
able bit rate from 6 to 54 Mbps. OFDM was also one of the promising candidates for the 
European third-generation personal communications system (universal mobile telecom¬ 
munication system). However, it was not approved since the code division multiple ac¬ 
cess (CDMA)-based proposals received more support. OFDM is now being considered 
for the fourth-generation mobile communication systems [3]. Therefore, OFDM’s per¬ 
formance in mobile and fading environments is the topic of many current studies. These 
studies are also helpful to determine the ways by which we can use this technique in bat¬ 
tlefield applications more efficiently. 

A. OBJECTIVE 

Although OFDM seems to be a solution to keep up with the demand of increasing 
data rates, it has some drawbacks. Sensitivity to frequency errors is the most significant 
of these drawbacks. The main objective of this thesis was to investigate and document the 
effects of frequency offset and phase noise on the performance of OFDM based digital 
communications under different channel conditions. 

A step-by-step approach was adopted in order to achieve the objective of this the¬ 
sis. The first step is to provide a basic background on the principles of OFDM. The rea¬ 
sons for the frequency errors and a theoretical analysis of these effects on OFDM systems 
are documented. Finally, a communication system utilizing OFDM was simulated in 
Matlab and the system’s behavior under different channel conditions was observed. 

To be able to observe the system behavior, the simulation results for different 
channel models are presented in graphical form. Next, the simulation results obtained in 
this work are compared to the simulation results reported in related studies. 
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B. RELATED RESEARCH 

Due to its many attractive features, OFDM has received much attention in the 
wireless communications research communities. Numerous studies have been performed 
to investigate its performance and applicability to many different environments. Below 
are some of the many studies conducted concerning the effect of frequency errors on 
OFDM systems. 

Many of the published studies about the frequency errors use two main refer¬ 
ences. The first is the study of Pollet el al. on sensitivity of OFDM systems to frequency 
offset and Wiener phase noise [4], and the second is the study of Moose on a technique 
for OFDM frequency offset correction [5]. 

Other related studies include the study of Armada on the phase noise and sub- 
carrier spacing effects on OFDM systems’ performance [6], the study of Xiong about the 
effect of Doppler frequency shift, frequency offset, and phase noise on OFDM receivers’ 
performance [7] and the study of Zhao el al. on the sensitivity of OFDM systems to Dop¬ 
pler shift and carrier frequency errors [8]. 

C. ORGANIZATION OF THE THESIS 

Chapter II starts with an introduction to OFDM and proceeds to present the fun¬ 
damentals of OFDM. The generation of an OFDM symbol is also explained using a step- 
by-step approach. Chapter El discusses the frequency errors including the frequency off¬ 
set and the phase noise. A simple OFDM signal representation is adopted to easily repre¬ 
sent the frequency offset and Doppler frequency shift effect in analytical expressions. 

The reasons for frequency errors are explained. 

Chapter IV develops the Matlab simulations to study the system performance un¬ 
der different environments. It also presents the results. Chapter V provides a summary of 
the thesis, the conclusions and suggestions for future work. 

Appendix A outlines the Matlab code for the transmitter, the channel and the re¬ 
ceiver. The function of each m-file and its relationship to other functions are discussed. 

At the end of each section, the hierarchy of the m-files is illustrated in a schematic form. 
Appendix B presents the Matlab code. 
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D. OFDM-BASED IEEE 802.11A STANDARD 


The rapid growth of the applications utilizing digital communication systems in¬ 
creased the need for high-speed data transmission. New multi-carrier modulation tech¬ 
niques are being proposed and implemented to keep up with the demand of higher data 
rates. Of these multi-carrier techniques, OFDM is the method of choice for high-speed 
communication due to its many attractive features. This chapter attempts to justify the 
choice of OFDM among other communication techniques. 

A. INTRODUCTION TO OFDM 

OFDM was introduced in the 1950s but was first implemented in the 1960s. It 
was originally developed from the multi-carrier modulation techniques used in high fre¬ 
quency military radios. A patent for OFDM was granted in the 1970s. However, when 
OFDM was first introduced, it was not very popular because of the cost and complexity 
of large arrays of sinusoidal generators and coherence demodulators [9]. 

The actual widespread use of OFDM started after the inverse discrete Fourier 
transform (IDFT) and discrete Fourier transform (IDFT) made the OFDM implementa¬ 
tion possible without the use of large number of sinusoidal generators. 

OFDM was accepted as a wireless local area network (WLAN) standard in Sep¬ 
tember 1999, but actually it was not the first IEEE physical standard for WLANs. The 
first standard was approved in June 1997 and specified one medium access control 
(MAC) and three physical layers (IEEE 802.11 FHSS, IEEE 802.11 DSSS and IEEE 
802.11 IR). The IEEE 802.11 direct sequence spread spectrum (DSSS) originally sup¬ 
ported both 1 Mbps and 2 Mbps of data rates while the other two supported 1 Mbps or 2 
Mbps optionally. Due to the increasing demand for higher throughput, a high-data-ratc 
DSSS (IEEE 802.1 lb) was selected for standardization in July 1998, and the data rate 
increased to 11 Mbps. The IEEE 802.1 la and IEEE 802.1 lb standards were developed 
simultaneously. The federal communications commission (FCC) released 300 MHz of 
spectmm in the 5.2-GHz band in January 1997 for WLAN applications, and IE EE 
802.1 la was targeted to use this spectral band [10]. 
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OFDM was chosen to be the physical standard for IE EE 802.1 la. The time- 
related parameters of OFDM for IE EE 802.1 la are listed in Table 1. 


Parameter 

Value 

Number of data sub-carriers 

48 

Number of pilot sub-carriers 

4 

Total number of sub-carriers 

52 

Sub-carrier frequency spacing 

0.3125 MHz 

IFFT/FFT period 

3.2 ps (1/ A f ) 

Preamble duration 

16ps 

Signal Duration BPSK-OFDM 
symbol 

4 Ms ( T gi +T fft ) 

Guard Interval (GI) duration 

0.8 ps (W 4) 

Training symbol GI duration 

1.6ps {T fft /2) 

Symbol interval 

4 Ms {T gi +T fft ) 

Short training sequence duration 

8ps (\0T fft /4) 

Long training sequence duration 

8 ms (T ci + 2T fft ) 


Table 1. OFDM Timing Related Parameters (From Ref. 11.) 


The rate dependent parameters of IEEE 802.1 la are listed in Table 2. As can be 
observed from Table 2, different sub-carrier modulation schemes can be used in combi¬ 
nation with different convolutional codes to achieve the desired data rate for the applica¬ 
tion. 


Data rate 
(Mbps) 

Modulation 

Coding 

Rate(R) 

Coded bits 
per sub¬ 
carrier 

(A t bpsc) 

Coded bits 
per OFDM 
symbol 

(/V CBPs) 

Data bits 
per OFDM 
symbol 

(A Ditrs) 

6 

BPSK 

1/2 

1 

48 

24 

9 

BPSK 

3/4 

1 

48 

36 

12 

QPSK 

1/2 

2 

96 

48 

18 

QPSK 

3/4 

2 

96 

72 

24 

16-QAM 

1/2 

4 

192 

96 

36 

16-QAM 

3/4 

4 

192 

144 

48 

64-QAM 

2/3 

6 

288 

192 

54 

64-QAM 

3/4 

6 

288 

216 


Table 2. OFDM Rate Dependent Parameters (From Ref. 11.) 
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Table 3 shows the major parameters of the OFDM physical layer. In the simula¬ 
tions conducted in this thesis, a data rate of 24 Mbps, modulation type of QPSK, and a 
coding rate of Vi were used, while the other parameters in Table 3 were left unchanged. 


Information data rate 

6,9,12,18,24,36,48 and 54 Mbps 
(6,12 and 24 Mbps are mandatory) 

Modulation 

BPSK OFDM 

QPSK OFDM 

16-QAM OFDM 

64-QAM OFDM 

Error correcting code 

K = 7 convolutional code 

Coding rate 

1/2, 2/3, 3/4 

Number of sub-carriers 

52 

OFDM symbol duration 

4ps 

Guard interval 

0.8 ps (T a ) 

Occupied bandwidth 

16.6 MHz 


Table 3. Major Parameters of OFDM Physical Layer (From Ref. 11.) 


OFDM is being used in a number of wired and wireless voice and data applica¬ 
tions due to its flexible system architecture. Some examples of current OFDM applica¬ 
tions are DAB (digital audio broadcasting), HDSL (High-Rate Digital Subscriber Line), 
VHDSL (Very High-Rate Digital Subscriber Line), ADSL (Asymmetric DSL), HDTV- 
Terrestrial (High Definition TV-T), IEEE 802.11 and HiperLAN/2, GSTN (General 
Switched Telephone Network), Cellular Radio and DVB-T (digital video broadcasting- 
terrestrial) [9]. 

From the range of the current applications and the research activity on OFDM and 
considering the fact that many wireless communication systems and network standards 
are being developed based on OFDM, it may be easily remarked that this technique will 
have a significant impact on the future of digital communications and wireless networks. 
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B. FUNDAMENTALS OF OFDM 


1. Multi-Path (Delay-Spread or Time Dispersion) 

In general, high data rate means short symbol time compared to the delay spread 

[Tsymbol < ^delay ) • Delay-spread greatly affects the communication system and the signal 
might not be recovered at the receiver. 

This section addresses the effects of delay spread which occurs as the surfaces be¬ 
tween a transmitter and a receiver reflect a transmitted signal as in Figure 1. 



Figure 1. A Multi-Path Scenario 


The receiver obtains the transmitted signals with random phase offsets and this 
causes random signal fades as reflected signals destructively or constructively affect each 
other [12], as seen in Figure 2. 


Transmitted signal 


r 
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Figure 2. Delayed Signals (After Ref. 13.) 
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When T symhol <T delay (B c <B S ), as in Figure 3, the signal faces frequency selective 

fading and this causes time dispersion. The effect of this is intersymbol interference (ISI), 
where the energy of one symbol leaks into another symbol, as can be viewed from Figure 
4. As a result, the bit error rate (BER) increases, this in turn degrades the performance. 

ISI is one of the biggest problems of digital communication and OFDM deals with this 
problem very effectively. 


Be 



(a) (b) 

Figure 3. Representation of a Symbol in A Frequency Selective Channel: (a) Time 

domain (b) Frequency domain 
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Figure 4. Illustration of ISI (After Ref. 14.) 
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A way to deal with frequency selective fading is to decrease the data rate and thus 
change the frequency selective fading to flat fading. The desired scheme is illustrated in 
Figure 5. OFDM systems mitigate the ISI by changing the frequency selective fading 
channel to flat fading channel as discussed below. 
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Figure 5. (a) Time Domain Representation, (b) Frequency Domain Representation of 

a Symbol in Flat Fading Channel. 

OFDM modulates user data onto tones by using either phase shift keying (PSK) 
or quadrature amplitude modulation (QAM). An OFDM system takes a high data rate 
stream, splits it into N parallel data streams and transmits them simultaneously. As can be 
observed from Figure 6, each of these parallel data streams has a rate of R/N , where R is 
the original data rate. The data streams are modulated by different carriers and combined 
together by inverse fast Fourier transform (IFFT) to generate the time-domain signal to 
be transmitted [12]. 
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Figure 6. OFDM Splits a Data Stream into N Parallel Data Streams (After Ref. 15.) 
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By creating a slower data stream, the symbol duration becomes larger than the 
channel’s impulse response. In this way, each carrier is subject to flat fading. 

2. Orthogonality 

Normally, multicarrier systems, such as frequency division multiplexing (FDM), 
have to modulate different sub-carriers with spectrally separate symbols to avoid inter 
carrier interference (ICI) at the cost of a bandwidth loss. However, in OFDM, spectrally 
overlapped sub-carriers can be used and since they are orthogonal, they do not interfere 
with each other. This makes OFDM a bandwidth efficient modulation scheme [12]. 

Orthogonality of the sub-carriers must be ensured to avoid ICI. OFDM carriers 
are orthogonal over a symbol interval only if they are spaced in frequency exactly at the 
inverse of the symbol period (see Figures 7 and 8). IFF ! implementation inherently pro¬ 
vides this frequency spacing requirement. It is also necessary to add a guard interval in 
multi-path channels to account for the multi-path spreading effects. Figure 7 shows a sub- 
carrier having a rectangular pulse shape and its frequency spectrum. The rectangular 
pulse has a width of T s and has a sine-spectrum with nulls at multiples of \/T s . 


Ko| 


Figure 7. Representation of an OFDM Sub-Carrier: (a) In Time Domain (B) Fre¬ 
quency Domain 

Figure 8 shows the spectra of multiple carriers. Although all frequency spectra are 
overlapping, they do not overlap at the carrier frequencies /,,/ 2 ,..., and this implies or¬ 
thogonality. In the time domain each carrier has an integer number of periods within the 
symbol time interval. Due to the orthogonality, the nulls of the side lobes coincide with 
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the peaks of the main-lobes of each carrier. By building the receiver to sample at the cen¬ 
ter frequency of each sub-carrier, it is possible to obtain only the corresponding signal’s 
energy, since the side-lobes of other carriers have zero energy at that frequency [12]. 



Figure 8. Sub-Carrier Orthogonality (After Ref. 16.) 


Orthogonality can be lost due to two main reasons. The first reason is the fre¬ 
quency errors, which are discussed in detail in Chapter in, and the second is the error in 
symbol boundary and sampling rate offset. 

3. OFDM Transmitter 

A block diagram of the OFDM transmitter module is presented in Figure 9. Each 
of the blocks is explained in detail in the following subsections. 
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Figure 9. Schematic Diagram of an OFDM Transmitter 
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a. Channel Coding 

A sequential binary input data stream is first encoded by the channel 
coder. Error correction coding is important for OFDM systems used for mobile commu¬ 
nications. When channel coding is used to improve its performance, OFDM is referred to 
as coded OFDM (COFDM). The IEEE 802.11a standard uses convolutional coding as 
forward error correction (FEC) coding, while the convolutional coder has rate Vi, con¬ 
straint length 7 and generator polynomials (133, 171) [11]. 

Table 2 lists the rate and modulation scheme combinations for a desired 
value of data rate. Higher data rates are obtained by puncturing the rate Vi code. The 
Convolutional encoder used in the IEEE 802.1 la standard is presented in Figure 10. 


Input 




► 

Output 

Figure 10. Convolutional Encoder with Constraint Length= 7. Each box represents a 

shift register (After Ref. 11.) 

b. Interleaving 

The encoded symbols, out of the convolutional encoder, are next inter¬ 
leaved by a block interleaver. The block interleaver takes data from the encoder, per- 
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mutes it and then outputs the rearranged binary data. The idea behind interleaving is to 
handle burst errors. The interleaver simply scatters the burst errors randomly so that the 
receiver can easily correct the random errors. 

Many interleaving methods exist, but the conventional block interleaver is 
widely used since it is the easiest to implement. To implement block interleaving, a ma¬ 
trix that can accommodate all the symbols in the transmit message matrix must first be 
created. Next, the symbols from the message matrix are read into this intermediate matrix 
by rows and the interleaving effect is created by reading these symbols out by columns. 
The degree of interleaving is determined by the dimensions of this intermediate matrix. 

c. Symbol Mapping 

For the symbol mapping part of the block diagram, the required data rate 
for the application will determine which signal constellation, i.e., which modulation tech¬ 
nique is to be used. 

Following the interleaver, the rearranged binary data are mapped onto the 
desired constellation points. Figure 11 shows two constellation examples that were listed 
in Table 2. 



Figure 11. 16-QAM and QPSK Constellations (From Ref. 11.) 
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d. Differential Quadrature Phase Shift Keying (DQPSK) Modula¬ 
tion 

After the symbol mapping process, the data can be frequency or time dif¬ 
ferential encoded. The Matlab implementation details of differential encoding are ex¬ 
plained in Appendix A. 


Differential (noncoherent) modulation techniques decode the incoming 
symbols by comparing the received symbol with the preceding symbol. If the transmitted 
waveform is given as [17] 


4 *) 



cos[vv+e,.(r)] 


i=l,...,M 0 <t<T, 


( 0 . 1 ) 


then the received signal is given by 


1 22J 

• (0 = j— cos [w 0 t + 0,. (t) +a ] + n(t) 


i = 0 <t<T 


( 0 . 2 ) 


where a is assumed to be a random variable uniformly distributed between 0 and 2% and 
n ( t ) is an additive white Gaussian noise (AWGN) process. 


For noncoherent detection, matched filters cannot be used because the 
matched filter output becomes a function of the unknown phase a . However, assuming 
that a changes slowly relative to the two symbol period times, the phase difference of 
two respective signals becomes independent of a and given by [17] 

[S,(7' 1 )+a]-[8 ; (r)+a]=6 t (7' 2 )-e i (r i ) = t,(r), (0.3) 

which suggests that the phase of the previous signaling interval can be used as a phase 
reference to demodulate the current signal. Of course, this requires the use of differential 
encoding in the transmitter. 

Even though DPSK performs less efficiently than PSK, it was chosen due 
to the ease of implementation. 

e. IF FT 

Hie interleaved and differentially encoded complex frequency array is 
next modulated onto a desired number of sub-carriers (52 sub-carriers are used in the 
TREE 802.1 la standard) by implementing an IFFT. 
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The OFDM baseband sub-carrier is 

§ k (0 = e j2nfk ' (0.4) 

where f k is the k ,h sub-carrier frequency. An OFDM symbol consists of iV modulated 
sub-carriers. The OFDM signal not including a cyclic prefix is given by [3] 

= 0 <t<NT (0.5) 

VA' k= 0 

where z , k is the k' h complex data symbol and NT is the OFDM symbol duration. The 
sub-carriers in Equations 2.4 and 2.5 have frequencies 


in the sense that 



( 0 . 6 ) 


, NT 

= (0.7) 

iVi 0 

ensures orthogonality. If the signal s(t) is sampled with a sampling period of T, the fol¬ 
lowing is obtained: 

s (n) = s(t = nT) = lDFT{z k } = -^Y,z k e J2nkn/N n = 0,l...N-l. (0.8) 

£ = o 

This last equation is IDFTjz, ] and was proposed by [18]. As can 

be seen from Equation 2.8, a baseband OFDM transmission symbol is an iV-point com¬ 
plex modulation sequence. It is composed of N complex sinusoids, which are modulated 
with z(k). 

f. Guard Interval Addition 

The addition of the Guard interval, T CI , is the next step in generating 
OFDM symbols. The guard interval is added to the start of OFDM symbols at the trans¬ 
mitter to e li minate ISI and ICI, and it is deleted before the EFT process at the receiver. 


Hie interval T C1 is an important parameter in the IEEE 802.1 la standard 
because it determines the choice of other parameters. It is normally desired that T C1 be 
larger than the expected multi-path delay spread to combat ISI. However, the larger T CI , 
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the smaller the effective symbol duration, thereby resulting in reduced throughput. In the 
IEE E 802.1 la standard, the total OFDM symbol duration is chosen to be 4 ps , and T CI is 
0.8 ps . Therefore, the effective symbol duration is 3.2 ps , and the inverse of the effec¬ 
tive symbol duration provides the sub-carrier spacing of 0.3125 MHz. These time related 
parameters of OFDM are listed in Table 1. 

The expression for a low-pass OFDM signal is given by Equation 2.8. The 
transmitted time domain signal is then obtained by including the guard interval. The ana¬ 
log signal after the digital to analog converter (DAC) in the block diagram is expressed as 

%«> = !<0 

A k=0 

where A/ is sub-carrier spacing. The guard interval Tgi is created by appending the IFFT 
sequence with a circular extension (called cyclic prefix) in order to maintain orthogonal¬ 
ity. For example, the cyclic prefix for the transmission sequence j s lp Or)} consists of 

the last m samples { s lp (N- m), s Ip (N - m -1),... s jp (N -1)}, which are added to the begin¬ 
ning of the sequence j .5^(0), s, (l),....?, (iV-1) j. The new sequence now consists of N+m 

samples. The benefit of the guard interval can be observed in Figure 12. Figure 12b illus¬ 
trates that the decaying transient of symbol 1 in the guard interval can be allowed since it 
does not affect symbol 2. 
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(a) OFDM Symbols without Guard Interval 
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(b) OFDM Symbols with Guard Interval 
Figure 12. Effects of Guard Interval (After Ref. 16.) 


The disadvantage of the guard interval is that it causes a loss of bandwidth 
efficiency since data are not transmitted during the guard interval. On the other hand, by 
adding the guard interval, ICI can be avoided since orthogonality is ensured, and ISI is 
avoided since adjacent symbols do not interfere with each other. 
g. Number of Sub-carriers 

The number of sub-carriers is another important parameter in the design of 
OFDM systems along with the system bandwidth B ~\/T . the sub-carrier bandwidth 

1 J NT , and the guard interval T CI . 

Initially, the guard interval T CI is chosen to be larger than delay spread x , 
and the OFDM symbol length should be larger than T CI . Therefore NT » x , where N is 

the number of sub-carriers which yields N » Bz . However, a large N means smaller fre¬ 
quency spacing between OFDM carriers, which in turn, makes the system vulnerable to 
Doppler spread and ICI. This last restriction forces larger frequency spacing than the 
Doppler spread (B/N » f d ) . As a result, the number of sub-carriers can be chosen in the 
interval given by [3] 
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B 

Tfi « iV « —. 

h 


( 0 . 10 ) 


The number of sub-carriers is fixed at 52 in the TERR 802.1 la standard. Pour of 
these sub-carriers are used as a reference to rrunirnize frequency and phase shifts of the 
signal during transmission. 

h. Symbol Pulse Shaping (Windowing) 

The time dispersion can be combated by the addition of guard interval. 
However, it is also necessary to combat the slow decay of OPDM symbols in the fre¬ 
quency domain. The effects of slow decay into the adjacent bands can be avoided by nar¬ 
rowing the spectrum. 

The method of choice to narrow the spectra is to use pulse shaping of the 
OPDM symbol. Pulse shaping is accomplished by using either a time window or a filter 
[3]. Windowing helps to smooth the transition between symbols and narrows the signal 
spectrum. 

i. RF Modulation 

After the pulse shaping process, the OFDM symbols are upconverted to a 
radio frequency carrier, amplified and transmitted through the antenna. The TPPP 
802.1 la standard specifies the 5-Ghz band for the carrier frequency. In this thesis, the 
communication system was simulated at the baseband. 

4. OFDM Receiver 

A block diagram of the OFDM Receiver module is presented in Figure 13. 
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Figure 13. Schematic Diagram of an OFDM Receiver 
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The processes performed in the transmitter are simply reversed in the receiver. In 
addition, the receiver must include algorithms for synchronization and channel estima¬ 
tion. In particular the synchronization block performs operations, such as packet detec¬ 
tion, timing estimation, frequency offset estimation and correction. 

a. Removing Guard Interval and FFT Processing 

Since only the baseband processing is being used in this thesis, the down- 
conversion process in the receiver is not needed. The first step is to remove the guard in¬ 
terval to obtain the information portion of the symbol for further processing. Next, the 
time domain samples are transformed into the frequency domain by the FFT process. 

This also makes it possible to recover the OFDM frequency tones. 

b. Decoding and Deinterleaving 

The next step in the receiver is the time or frequency differential decoding. 
Following the differential decoding, the inverse mapping of each received complex 
modulation value into a corresponding N- ary symbol is accomplished. The message is 
next deinterleaved according to the particular interleaving configured in the transmitter. 

c. Received message and Viterbi Decoding 

After being exposed to the effects of different channel models, some errors 
may occur in the transmitted message. The received bit stream is determined after the 
Viterbi decoder. 

The Viterbi decoders are used to decode bit streams encoded by convolu¬ 
tional encoders. After receiving the bit stream from the channel, the Viterbi algorithm 
searches the tre ll is to find the path that can create the bit stream closest to the received 
one. The measure of closeness is different for hard and soft decision decoding. For the 
hard decision decoding, the Viterbi algorithm tries to find the bit stream path that has the 
minimum Hamming distance from the received sequence. For soft decision decoding, the 
algorithm tries to find the bit stream path that has the minimum Euclidian distance from 
the received sequence. The details of the Viterbi decoding algorithm can be found in [17]. 
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C. SUMMARY 

An introduction to OFDM was presented in this chapter. The transmitter and re¬ 
ceiver block diagrams have been discussed in detail and the chapter concluded by listing 
the advantages and disadvantages of OFDM and mentioning the potential problematic 
areas in OFDM implementation. 

The advantages of OFDM can be summarized as follows. OFDM changes fre¬ 
quency selective fading to flat fading and offers high bit rate transmission. Its flexible 
architecture allows the use of combinations of various coding schemes (rates) and modu¬ 
lation schemes (signal constellations) according to channel condition. OFDM e lim inates 
ISI, is bandwidth efficient, and supports several multiple access schemes (TDMA, 

FDMA, MC-CDMA, etc.) [15]. 

On the other hand OFDM is very sensitive to frequency errors, has a large peak to 
average ratio (PAR), and is sensitive to channel fading [15]. The effect of PAR is ana¬ 
lyzed in [19]. 

OFDM has some potential problematic areas that should be considered during 
implementation. The first is synchronization, including the time and frequency synchro¬ 
nization. The second is the non-constant power envelope, which requires the use of linear 
amplifiers. The last problematic area is the requirement of channel estimation since the 
channel is time-variant [13]. 

The next chapter presents the steps to create an OFDM signal. OFDM systems’ 
unique requirements, such as IFFT/FF1 algorithms and the cyclic prefix, are also cov¬ 
ered. Following that, the effects of frequency errors are analyzed. 
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ID. FREQUENCY ERRORS IN OFDM AND THEIR EFFECTS 


Before an OFDM receiver demodulates the sub-carriers, it performs two synchro¬ 
nization tasks. First, it determines the symbol boundaries and the optimal timing instants 
in order to min im ize both the ICI and the ISI. Second, it tries to estimate and correct the 
frequency errors [20]. OFDM’s sensitivity to these frequency errors is one of its main 
drawbacks. This chapter discusses the effects of these frequency errors on OFDM sys¬ 
tems. 

A. FREQUENCY ERRORS 

The orthogonality of the sub-carriers can be ensured only if the receiver and the 
transmitter have the same reference frequency. Any deviation from this reference fre¬ 
quency may cause ICI and loss of orthogonality. Another frequency error factor is phase 
noise, which is caused by random jitter of the phase of the steady sinusoidal waveform 
generated by the oscillators [20]. 

Typically frequency errors are generated by the fact that the oscillators in the 
modulator and demodulator do not have exactly the same frequency. 

For single-carrier systems, the effect of phase noise and frequency offset appear 
only as degradation in the received SNR, rather than ISI or ICI. Nevertheless, many effi¬ 
cient techniques to minimize the effects of this drawback have been proposed in the lit¬ 
erature [20, 21]. 

Other reasons for frequency errors include Doppler shift caused by the relative 
movement between the receiver and the transmitter, and phase noise introduced by non¬ 
linear channels. 

Figure 15 shows the front end of an OFDM receiver where most of the frequency 
errors occur, i.e., the local oscillators and the sample clock at the analog to digital (A/D) 
converter [3]. The A/D converter causes errors when the receiver does not have the same 
sample clock frequency as at the transmitter. 


21 




(LO) f 

sample clock 

Figure 14. OFDM Receiver Front End (After Ref. 3.) 


The frequency offset and the phase noise cause a phase rotation of the received 
symbols. Coherent OFDM systems need a phase tracking device to obtain the phase of 
the incoming symbols for correct demodulation [22]. Typically, three types of algorithms 
are used to estimate the frequency offsets, i.e., to track the phase, in coherent OFDM sys¬ 
tems: 

• data-aided algorithms that use special training information (pilot tones) 
within the transmitted signal, 

• non-data-aided algorithms that analyze the received signal in the fre¬ 
quency domain, and 

• algorithms that use the cyclic prefix of the OFDM signals [21]. 

For wireless applications, the first method is the method of choice. The imple¬ 
mentation details of these coherent OFDM detection algorithms can be found in [21]. It 
should be noted that in many applications as well as in this thesis, the differential detec¬ 
tion technique is used since it does not require channel estimation and the use of pilot 
tones, even though it causes a loss of SNR. 

In the following sections, the effects of each of the aforementioned reasons that 
cause frequency errors are analyzed. A simple representation of an OFDM signal is cre¬ 
ated for the ease of analysis. 
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B. CREATING A SIMPLE REPRESENTATION OF AN OFDM SIGNAL TO 
SIMPLIFY ANALYSIS 

A baseband OFDM signal can be represented by [7] 

N -1 

b(t)=Y,A cos ( ( ®i t+ §i) ( 0 . 11 ) 

;=1 

where A i is the amplitude, co, = 2k J] is the angular frequency, <f» i . is the phase of the i ,h 
sub-carrier, and N is the number of sub-carriers. 


According to the modulation technique to be used, either A or (f) is determined by 
the data. As discussed in Chapter n, the baseband OFDM signal b(t) is modulated next, 
onto a RF carrier with frequency f c : 


s(t ) = 2b(t)cos(O c t 

N -1 

= 2^2 A c os(( 0 ,t + A)cos (a c t 

i=0 

N -1 

= E A { cos [( ro c +©/ )t+h] +COS [(<D C -CO, ) t - A]} 

i=0 


( 0 . 12 ) 


where ot =2nf c , and we assume the phase of the carrier to be zero for simplicity. Since 
a single side band transmission is enough to carry the information in A i or (]),, it is as¬ 
sumed that the upper sideband is used, and therefore the transmitted signal can be repre¬ 
sented as [7] 

JV—1 

s(t) = 52 A cos[(co c +( 0 , )f+()),]. (0.13) 

i=0 

C. THE EFFECT OF DOPPLER FREQUENCY SHIFT, PHASE NOISE AND 
OSCILLATOR FREQUENCY OFFSET ON THE OFDM SYSTEMS 

In this section the theoretical analysis of the effects of frequency errors is pre¬ 
sented. The maximum Doppler shift occurs when the two mobile nodes move toward 
each other, given by [23] 
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(0.14) 



c 

where v is the relative speed of the two nodes, f c is the carrier frequency and c is the 


speed of light(3xl0 8 m/s). 

1. Effect of Doppler Shift on the Carrier Frequencies, Sub-carriers, En¬ 
velope and Symbol Timing 

An OFDM signal consists of numerous sub-carriers with different frequencies. 

The amount of Doppler shift affecting the i ,h sub-carrier is given by [7] 

DopP ' er slufl > (1 + ^){f c ±f) (0.15) 

where c, is the percentage of the change in frequency and is determined by 

£, = —= —cosO. (0.16) 

/ c 

Tlie right-hand side of Equation 3.5 can be written as 

(i+4)(/c±/,)=(i+5)/c±(i+5)/, (0-17) 

which demonstrates that the Doppler frequency shift affects the carrier frequency and the 
sub-carrier frequencies by the same percentage c, [7]. The Doppler shift of the carrier 

frequency can be calculated as 

fdc = —cose (0.18) 

c 

and the Doppler shift of the sub-carrier frequencies as 

fa=— cosO- (0.19) 

c 

By using Equation 3.3 again, the transmitted OFDM signal with Doppler shift can 
be written as 


N -1 


s(t) = £a, cos[( 1+£)(cq +co, )t+f] 


/=o 

N -1 


( 0 . 20 ) 


= E { A c °s [( 1 + 4) 0)/ + 4>,. ] COS [ (1 + %) 0)/] - A,, sin [( 1 ^ + £) co/ + <> f ] sin [( 1+^) m /]} . 

i =0 
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In Equation 3.10, A. cos [(1 + q )coi + (J) ; ] can be thought of as the envelope of the 

carrier, cos [(1 + ^) co ?] , which helps to demonstrate that the Doppler shift affects the 

envelope and the carrier frequency by the same percentage [7]. The Doppler shift also 
affects the symbol rate and the time synchronization. 

2. Phase Noise, Oscillator Frequency Offset and Their Effects 
The carrier frequency offset has two destructive effects on OFDM systems. The 
first is the reduction of the amplitude of the desired sub-carrier and the second is the in¬ 
troduction of ICI [21]. The reduction of the amplitude happens because the desired sub- 
carrier is not sampled at the peak of the sine function of the DFT since the sine functions 
are shifted (see Figure 15). Also, the ICI is caused since orthogonality is lost between the 
neighboring sub-carriers [24]. 


Frequency 



Moose [5] and Pollet el al. [4] analytically evaluated the effects of the carrier fre¬ 
quency offset and carrier phase noise on the SNR degradation for an AWGN channel. 
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The results derived in these papers are used for many of the studies on frequency errors. 

As a result, we have used these results as benchmar ks for the simulation results in this 
thesis. 

Moose [5] models the channel by a complex transfer function H k for the k' b sub- 
carrier and assumes that a frequency offset of £ exists for all sub-carriers. The frequency 
offset, £ , is usually expressed relative to the data symbol rate, and it includes both the 
Doppler shift and the local oscillator frequency offset effects. As mentioned previously, 
each sub-carrier is affected by the Doppler shift according to its frequency value. Taking 
all the causes into consideration, the total offset for each sub-carrier after down conver¬ 
sion can be expressed as A f t = f lc + f di + Af [0 where f dc is the carrier Doppler shift, f di is 

the i' h sub-carrier Doppler shift and A f L0 is the local oscillator frequency offset [7]. 


As can be seen from Equations 3.8 and 3.9, f dc » f di since f c » f ]. Thus, 

A f t = f dc + A f L0 = A/ . This result demonstrates that the frequency offset is independent of 
the sub-carriers, which in turn states that the relative offset £ = A f /R s is also independ¬ 
ent of sub-carriers [7]. 

The SNR calculated at the output of the DFT of the receiver is given by [5] 


SNR > 


N„ 


sin TIE 
7t£ 


1 


J 1 + 0.5947 




(sin7t£ )' 


|e|<0.5 


where EjN 0 is the OFDM carrier energy to the AWGN spectral density. An upper 
bound on the degradation can be obtained from Equation 3.11 as follows: 


Ab -lOlog 


10 


f E , 3 

1+ 0.5947 -^sin 2 7tA/ 
No 

sin c 2 Af 


( 0 . 21 ) 


( 0 . 22 ) 


V ) 

where the factor 0.5947 is derived by taking the lower bound of the summation of all in¬ 
terfering sub-carriers [24]. 
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Similar results to Equations 3.11 and 3.12 are found in [4] in which the channel is 
modeled by a time-varying phase 0 (t ), which is a result of either the phase noise of the 
carriers or a carrier offset between the receiver and transmitter. A distinction is made be¬ 
tween phase noise and frequency offset. 

a. Phase Noise 

For the phase noise case, 0 (t) is assumed to be a Wiener process with 
£"{0 (/)} = 0 and£j(0 [t o + t)-Q [t o ))‘j = 47t[5 |f|, where (Hz) is the one-sided 3-dB 
linewidth of the Lorentzian power density spectrum of the free-running carrier generator 

[4]. 


Phase noise has two main effects. First, it causes a random phase variation 
common to all sub-carriers. The effects of this common phase error are minimized by 
employing phase tracking techniques or differential decoding. Second, it introduces ICI. 
Based on the model defined in [4], the degradation D in SNR, i.e., the required increase 
in SNR to compensate for the phase noise is 


A® 


1 1 


^ An N — ^^' 
R 


,N 

J o 


6-lnl0| 

Since R= N/T = NR S , where N is the total number of sub-carriers and R s is the sub- 
carrier symbol rate, Equation 3.13 can be rewritten as 


(0.23) 


D m = 


11 


dB 


471 - 5 - 

R 


N n 


6-lnlO 

v J y “ 

The plot of SNR degradation in dB as a function of the norma liz ed band¬ 
width (p / /? ) is shown in Figure 18 for QPSK. This figure is plotted for 

EjN 0 = 10.5 dB corresponding to a BER of 10 6 for uncoded QPSK [20]. 


(0.24) 
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-3 dB linewidth/subcarrier spacing ^g-3 

Figure 16. SNR Degradation Versus - 3-dB Bandwidth of the Phase Noise Spectrum 

for QPSK. 

By observing Figure 16, it can be concluded that for a negligible SNR 
degradation of about 0.1 dB, the - 3-dB phase noise bandwidth should be approximately 
0.1 percent of the sub-carrier spacing for QPSK. However, this value changes depending 
on the modulation. 

b. Frequency Offset 

For the frequency offset case, 0 (t) is given by 

0 (t) = 2nAFt+Q 0 (0.25) 

where A F is the carrier offset. If there is a frequency offset, then the number of cycles in 
the FFT interval is not an integer and this causes ICI after the FFT. The sub-carriers in 
the middle of the OFDM spectrum are more exposed to ICI than the sub-carriers on the 
edges of the spectrum because they have interfering sub-carriers on their both sides [20]. 
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Based on the model developed in [4], the degradation D in SNR for fre¬ 
quency offset is given by 

10 


D ao = 


dB 


3-lnl0 
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(0.26) 


By using the relation R = N/T = NR S , Equation 3.16 can be rewritten as 
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where A F/R s is the relative frequency offset. The plot of SNR degradation in dB as a 
function of the relative frequency offset (A F/R s ) is shown in Figure 17 for QPSK. This 
figure is plotted for E s /N 0 = 10.5 dB corresponding to a BER of 10 6 for uncoded QPSK 
[ 20 ], 


(0.27) 



Figure 17. SNR Degradation Versus the Normalized Frequency Offset for QPSK. 
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From Figure 17, it can be seen that for a negligible SNR degradation of 
about 0.1 dB, the tolerable frequency offset should be approximately 2.5% of the sub- 
carrier spacing for QPSK. However, this value changes depending on the modulation. 

D. SUMMARY 

This chapter discussed the frequency errors and their causes. It was demonstrated 
that SNR degradation due to the Doppler shift, most of the time, is negligible, and this 
left frequency offset and phase noise of local oscillators as the main sources of degrada¬ 
tion. It was also shown that this degradation varies strongly according to the modulation 
used. In [21], it was illustrated that 64-QAM could not tolerate more than 1% carrier fre¬ 
quency error for a SNR degradation of 0.5 dB, whereas QPSK modulation could tolerate 
5% carrier frequency error for the same SNR degradation. It can also be stated that close 
spacing of carriers in frequency made the tolerable frequency offset a small percentage of 
the channel bandwidth. In [5], it is illustrated by simulation that the carrier frequency off¬ 
set is lim ited to 4% or less of the intercarrier spacing to maintain signal-to-interference 
ratios of 20 dB or greater for the OFDM carriers. 

In the next chapter, this thesis will attempt to simulate an OFDM system under 
frequency offset and phase noise conditions to observe the expected performance degra¬ 
dation discussed in this chapter. 
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IV. OFDM SYSTEM SIMULATIONS IN MATLAB 


The previous chapters covered the principles of OFDM and discussed the fre¬ 
quency errors that might occur in a communication system as well as the effects of these 
frequency errors. This chapter presents simulation of an OFDM communication system, 
operating under different channel conditions. The results of simulation are compared to 
the theoretical results discussed in Chapter IH. The simulation steps performed are pre¬ 
sented in Figure 18. 



Figure 18. Matlab Simulations Scheme 

A. SELECTING SYSTEM PARAMETERS 

Hie simulations attempt to emulate the bit rate and performance requirements of 
the IEEE 802.11a standard by using the Matlab software. An operational bit rate of 24 
Mbps was chosen and, according to the standard, the occupied bandwidth was taken as 
16.6 MHz. The constellation to be used was decided based upon these parameters. First, 
the number of bits/Hz needed to be transmitted to meet the requirements is calculated as 

— = 24MbpS = 1,44bit/Hz « 2bit/Hz. 

Hz 16.6MHz 
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QPSK signaling was chosen for this simulation, and the actual bit rate was 

2x16.6 = 33.2Mbps. Even though this seems much more than the required bit rate, the 

addition of the guard interval and EEC reduces the bit rate to the desired level [25]. 

The required symbol interval in the standard was given as 4 ps , and the sub- 
canier spacing was given in Table 1 as 0.3125MHz . The number of sub-carriers used in 

the standard is calculated using this sub-carrier spacing and occupied bandwidth as fol¬ 
lows: 

Ar occupied BW 16.6MHz 

N- ---=-= 53. 

frequencyspacing 0.3125MHz 

The standard actually uses 52 sub-carriers with 4 of these reserved for pilot sig¬ 
nals in coherent detection systems to make the system robust against frequency errors. 
Since differential encoding and decoding were used in this work, all 52 sub-carriers could 
be used to carry information symbols. 

A 16-point guard interval was used during the simulations to represent an 800-ns 
period, and a 64-point FFT was used to represent the 3.2 ps information signal length, 
i.e., a total of 80 time samples representing the standard 4 ps OFDM symbol. 

B. SIMULATION METHODOLOGY 

After modifying the code from [25, 26] and adding required sub-blocks to imple¬ 
ment the OFDM system for this thesis, a number of system simulations were performed 
using different channel models. 

The simulation process was performed in five steps. The first step was to verify if 
the Matlab code was properly functioning. After verification of the code, a channel with 
only AWGN was simulated, and the system’s performance in this channel was investi¬ 
gated. The third step was to investigate the performance of four mobile channels. During 
this step, how different mobile channels affect the system was demonstrated and, then, 
the most severe of them was chosen as the mobile channel. The first three simulation 
steps were conducted to ascertain the system’s response to different channel models 
without introducing any frequency offset and phase noise. The main goals of the simula¬ 
tions were achieved in the fourth and the fifth steps. The fourth step included the fre- 
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quency offset, which was studied for three different channels. In the fifth step, the phase 
noise was introduced and the performance measured for three different channels. Table 4 
summarizes all of the simulation steps. 


Step No. 

Simulation 

Channel Description 

1 . 

Code check 

Channel without AWGN and 

multi-path effects (Ideal channel). 

2. 

Performance in AWGN 

AWGN channel. 

3. 

Performance in Mobile 

channel 

Multi-path+AWGN channel. 

4. 

Performance with Fre¬ 
quency offset 

a. Ideal channel. 

b. AWGN channel 

c. Mobile channel 

5. 

Performance with Phase 

noise 

a. Ideal channel. 

b. AWGN channel 

c. Mobile channel 


Table 4. Simulation Steps 


This step-by-step simulation methodology, beginning with the individual channel 
simulations and advancing to the combined channel simulations helped to evaluate each 
channel model’s output individually. This approach made it possible to detect the errors 
caused by incorrectly designed sub-blocks at an early stage and forced a redesign of the 
incorrectly functioning blocks. 

The implementation details and results of each simulation step are explained be¬ 
low. 
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1. Code Check (Channel Model 0) 

The first step in the simulation process was to check whether or not the system 
model functioned satisfactorily in a noise-free, ideal channel. Due to the lack of any 
channel distortions affecting the signal, the transmitted and received messages should be 
exactly the same in the case of a correctly implemented overall system. After many simu¬ 
lations, with different parameters each time, it was determined that the system was func¬ 
tioning correctly since it was possible to receive the transmitted signals. 

Table 5 lists the parameters of a simulation run for QPSK realization. The trans¬ 
mitted and received QPSK signal constellations under noise-free channel conditions are 
shown in Figure 19, which verified a correctly functioning system since we were able to 
receive exactly what we had transmitted. 


# of car¬ 
riers 

# of sym¬ 
bols 

Channel model 

seed 

Interleaver ma¬ 
trix 

48 

300 

0 

22 

[28 12] 


Table 5. Code Check Simulation Parameters 
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Figure 19. (a) Transmitted and (b) Received QPSK Constellations After Differential 

Decoding, from a Channel Model 0 Simulation 
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2. OFDM Performance in an AWGN Channel (Channel Model 1) 

In this step, AWGN channel simulations were conducted, and the system’s per¬ 
formance results were compared with the reported values in other related studies [25, 26 
and 27]. 


The theoretical symbol error rate P b for BPSK is given by 


P b = Q 



(0.30) 


where [17] 


G(x,= ^f exp 


v' 2 / 


du. 


For M-ary modulation, P b can be approximated by assuming that each symbol error 
causes only one bit error as follows: 


P P 

s _ s 


log? M 


(0.31) 


where k is the number of bits and M = 2 k is the size of the constellation. For QPSK the 
symbol error rate is given by [21] 


P s = 2 Q 
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1 4e 
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(0.32) 


For higher order PSK modulation P s can be approximated by 


P s =2Q 



N 

V » 0 J 


sin 


( n 3 

k Mj 


(0.33) 


where E = E h (log? M ) is the energy per symbol. As seen from Equation 4.4, the bit er¬ 
ror rate can always be derived by dividing P s by the number of bits. 


An OFDM system using QPSK and the parameters listed in Table 6 was simu¬ 
lated under AWGN conditions. The received QPSK constellation is shown in Figure 20 
for a a value of 0.02. 
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# of car¬ 
riers 

# of sym¬ 
bols 

Channel 

model 

seed 

Interlever 

matrix 

<3 range 

48 

20000 

1 

33 

[144 139] 

0-0.08 


Table 6. AWGN Channel Simulation Parameters 
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Figure 20. Effect of AWGN on Received Signal Constellation for a = 0.02 (Before 

Differential Decoding) 

Figure 21 shows the simulation BER plot under AWGN conditions. Clearly, the 
errors decrease as the E b /N 0 increased. This result is used in later simulations as a refer¬ 
ence for comparing the results with frequency errors (see Figures 29, 30, 34 and 35, dis¬ 
cussed later in the thesis). 



36 










































Figure 21. BER versus E b /N a Plot for a QPSK Simulation in AWGN Channel 

3. OFDM Performance in Mobile Channel (Channel Model 2) 
a. Mobile Channels 

Mobile channels are composed of both multipath and AWGN components. 

The multipath component of mobile Channels 1 and 2 has a delay line with 8 taps. The 
input time domain signal is filtered by 4 parallel 7 th order FIR filters, whose coefficients 
are fixed according to [28]. These filters implemented in m-file cvdd as in Figure 22. The 
FIR filter output is computed using 

y i { n ) = ^ J a jkx{n-k) i = 1,2,3,4 (0.34) 

k =0 

where a ik are the filter coefficients and x(n) are the time-domain input samples. The 
filter coefficients used for each of the four FIR filters are fisted in Table 7. 
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k 

0 

1 

2 

3 

4 

5 

6 

7 

«i k 

- 0.013824 

0.054062 

- 0.157959 

0.616394 

0.616394 

- 0.157959 

0.054062 

0.013824 

a 2k 

0.003143 

- 0.019287 

0.1008 

- 1.226364 

1.226364 

- 0.1008 

0.019287 

- 0.003143 

a 3k 

0.055298 

- 0.216248 

0.631836 

- 0.465576 

- 0.465576 

0.631836 

- 0.216248 

0.055298 

«4 k 

- 0.012573 

0.077148 

- 0.403198 

0.905457 

- 0.905457 

0.403198 

- 0.077148 

0.012573 


Table 7. FIR Filter Coefficients for Mobile Channels 1 and 2 

■r(n) 



Figure 22. Modeling of Mobile Channels 1 and 2 Using FIR Filters and Delay Ele¬ 
ments D. 

Multipath components of mobile channels-3 and 4 use 6-tap FIR filters but 
they have time-varying coefficients with Jakes spectrum [29]. The filtered output is de¬ 
rived simply by summing the delayed paths with the direct path. 

All multipath components also include power loss of received signal levels 
and Doppler frequency shifting characteristics. (The delay, power loss and Doppler vec¬ 
tors can be altered to reflect the desired channel characteristics.) 

b. Mobile Channel Simulations: 

The goal of this simulation step was to see the mobile channel effects on 
the transmitted symbols. The four different mobile channel models were simulated and 
the performance curves for each of them presented. The one that causes the most degra¬ 
dation was chosen to be able to simulate the characteristics of a mobile outdoor environ¬ 


ment. 
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Mobile channel simulations were conducted with the parameters in Table 
8. Figure 23 shows the received signal constellations after each mobile channel before 
differential decoding for a noise variance value of a = 0.01 . Since the mobile channel 
causes constructive and destructive effects on the OFDM symbols, the received constella¬ 
tion points are scattered from their transmitted positions. If the differential encod¬ 
ing/decoding were not implemented, a large number of received symbols would be de¬ 
coded in error since many points are scattered into adjacent phase sectors. 


# of carri¬ 
ers 

# of sym¬ 
bols 

Channel 

model 

seed 

Interlever 

matrix 

<3 

48 

20000 

21,22,23,24 

33 

[144 139] 

0.01 


Table 8. Mobile Channel, Simulation 1 Parameters 
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(c) (d) 

Figure 23. Effect of Mobile Channels on Received Signal Constellations (Before Dif¬ 
ferential Decoding), (a) Mobile Channel 1 (b) Mobile Channel 2 (c) Mobile 
Channel 3 (d) Mobile Channel 4. 

Figure 24 shows the received QPSK constellations after differential de¬ 
coding. The scattered points are placed within their previous phase sectors, thereby re¬ 
ducing the number of errors. Figure 25 illustrates the effect of the mobile channels on the 
magnitude plot of the received symbols. As previously discussed, the sub-carriers in the 
middle of the OFDM spectrum experience more multipath loss than the sub-carriers on 
the edges of the spectmm because the sub-carriers in the middle have interfering sub¬ 
carriers on both sides. The effects of multipath can be observed in the form of magnitude 
distortion in Figure 25. Also, we can see the randomly generated peaks, which can be at¬ 
tributed to the use of the Jakes spectmm for mobile Channels 3 and 4. 
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(C) (d) 

Figure 24. Received Signal Constellation Affected by Mobile Channels (After Differ¬ 
ential Decoding), (a) Mobile Channel 1 (b) Mobile Channel 2 (c) Mobile 
Channel 3 (d) Mobile Channel 4. 
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(c) (d) 

Figure 25. Effect of Mobile Channels on the Magnitude Plot of Received Symbols, (a) 
Mobile Channel 1 (b) Mobile Channel 2 (c) Mobile Channel 3 (d) Mobile 

Channel 4. 

BER plots of the OFDM system under the mobile channel conditions were 
obtained by configuring the system with the parameters in Table 7. A range of a values 
from 0 to 0.08 was used to obtain these plots. 

Figure 26 shows the BER versus E b /N o plot of each mobile Channel. By 

observing the BER curves, it can be seen that mobile Channel 1 and 2 required more sig¬ 
nal energy relative to the other channel models for the same BER. Mobile Channels 1 and 
2 had approximately the same performance, and we randomly chose mobile Channel 1 as 
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the mobile channel model. This figure justifies that Doppler frequency shift have negligi¬ 
ble effects on the performance of OFDM systems since we get similar BER curves for 
mobile Channel 1 and 2. The Doppler frequency for mobile Channel 1 was 15 Hz and for 
mobile Channel 2 it was 5 Hz, according to [29]. 



Figure 26. BER versus E b /N a Plots for QPSK Simulation in Mobile Channels. 

4. Performance with Frequency offset 
a. Noise Free Channel 

This step of the simulation process was to verify if it were possible to im¬ 
plement the frequency offset. As previously mentioned in Chapter III, for the frequency 
offset case, 9 (Y) is deterministic and is calculated from27t A F t +0 O , where A F is the car¬ 
rier offset. In the literature, however, frequency offset is usually expressed as a ratio of 
the frequency spacing (AF/R s ). The performance degradation caused by frequency oflset 
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is also expressed in terms of this relative frequency offset. Therefore, in the simulations, 
the same approach was used, and each time domain sample of the OFDM symbols was 
exposed to the effective value of the relative frequency offset. 

To see the frequency offset effects exclusively, a noise-free channel was 
needed. Thus, Channel 0 was chosen. The easiest way to observe the effect of frequency 
offset was to compare the signal constellations of the transmitted and received symbols. 
The code to simulate a channel with only frequency offset was configured with the pa¬ 
rameters in Table 9. A relatively small frequency offset value was chosen in order to ob¬ 
serve its effect clearly. 


# of car¬ 
riers 

# of sym¬ 
bols 

Channel 

model 

seed 

Interleaver 

matrix 

A FIR, 
(%) 

48 

8000 

0 

22 

[167 48] 

0.1 


Table 9. Frequency Offset in Noise Free Channel Simulation 


Hie received QPSK signal constellation before differential decoding is 
shown in Figure 27 with the frequency offset. As seen in the figure, the center of the con¬ 
stellation points is shifted from their original positions. Figure 28 shows that the differen¬ 
tial decoding algorithm successfully placed the scattered points within their previous sec¬ 
tors. Due to the small frequency offset value, the errors were easily corrected. However, 
for wide ranges of frequency offsets, as the frequency offset value increased, so did the 
errors. 
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Figure 27. Received Signal Constellation in Channel Model 0 with Frequency Offset of 
0.1% of Frequency Spacing (Before Differential Decoding) 



Re 

Figure 28. Received Signal Constellation in Channel Model 0 with Frequency Offset of 
0.1% of Frequency Spacing (After Differential Decoding) 
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b. AWGN Channel 

This simulation step is important because all theoretical expressions for 
frequency offset, in the literature, were derived under AWGN channel conditions. Thus, 
the goal was to study how the system’s performance deviated from the theoretical per¬ 
formance with increasing frequency offset. 

The OFDM system was configured with the parameters given in Table 10. 


# of car¬ 
riers 

# of sym¬ 
bols 

Channel 

model 

seed 

Interleaver 

matrix 

G 

range 

AF/R S 

(%) 

48 

20000 

1 

33 

[144 139] 

0-0.8 

0-0.4 


Table 10. Frequency Offset in AWGN Channel Simulation 


Figure 29 show the BER performance for relative frequency offset values 
ranging from 0 to 0.4% in 0.1% increments. As can be seen from the figure, above 0.4% 

frequency offset BER reaches 10 1 , which is unacceptably high. 



Figure 29. BER versus E h /N o Plots for QPSK in AWGN Channel with Relative Fre¬ 
quency Offset Values from 0% to 0.4% 
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In Figure 29, the required E b /N a values increased paraboloically as the 
relative frequency offset values were increased; we compare these with the BER plot for 
AWGN simulation. This was an expected result because, by looking at Equation 3.17, the 
SNR degradation increased with the square of the relative frequency offset. 

c. Mobile Channel 

The mobile channel models were intentionally chosen to be severe chan¬ 
nels to represent the outdoor mobile environments. In this step, the goal was to see the 
frequency offset effect in a severe mobile channel. Channel 1 is selected for this reason. 
The system was configured with the same parameters presented in Table 9 but mobile 
Channel 1 was used for this simulation. 

Figure 30 shows the BER plot of the system in mobile channel 1 for rela¬ 
tive frequency offset values ranging from 0 to 0.4%. 



Figure 30. BER versus E b /N a Plots for QPSK in Mobile Channel 1 with Relative Fre¬ 
quency Offset Values from 0% to 0.4% 
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In Figure 30 the dotted lines are from a frequency offset simulation in the 
AWGN channel and the solid lines are from a frequency offset simulation in mobile 
channel 1. We can observe the degradation caused by the mobile channel by comparing 
the dotted and the solid lines and also the OFDM system’s performance in mobile Chan¬ 
nel 1 can be seen. 

5. Performance with Phase Noise 
a. Noise-Free Channel 

As the first step, the noise-free channel was used to study the phase noise 
effect exclusively. As mentioned previously, for the phase noise case, 0 (t) was assumed 

to be a Wiener process with ii{0 (r)] = 0 and £ j(0 (/ ; + /) -0 [t o )) j - 471 [11/|, where 

|j (Hz) is the one-sided 3-dB line width of the Lorentzian power density spectrum of the 
free-running carrier generator. A Wiener process was created in Matlab with zero mean 
and variance 4n\i |r| with |j as a percentage of the frequency spacing. The first simula¬ 
tion was run with the configuration parameters given in Table 11. Figure 31 shows the 
generated phase noise as a Wiener Process. 


#of 

carriers 

# of sym¬ 
bols 

Channel 

model 

seed 

Interlever 

matrix 

P/R, 

<%) 

48 

8000 

0 

222 

T167 481 

0.01 


Table 11. Phase Noise in Noise Free Channel Simulation 
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Figure 31. Phase Noise Generated as a Wiener Process 

Two different phase-noise sub-blocks were created. The first began the 
noise generation process with the first OFDM symbol and ended when the last time do¬ 
main sample of the last OFDM symbol transmitted, and the second created phase noise 
for each OFDM symbol individually, i.e., one Wiener process for each 80 time domain 
samples. The first method was used during the simulations since it more accurately re¬ 
flected the practical operation of a communication system. 

Figure 32 shows the received QPSK signal constellation as a result of 
phase noise, before differential decoding. Figure 33 shows the corrected signal constella¬ 
tion after differential decoding. The differential decoding algorithm successfully placed 
the scattered points within their respective sectors. 
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Figure 32. Received Signal Constellation in Channel Model 0 with fj / R = 0.0001. 

(Before Differential Decoding) 



Figure 33. Received Signal Constellation in Channel Model 0 with fj / R = 0.0001. 

(After Differential Decoding) 
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b. AWGN Channel 

To study the phase noise effect on the system, the code was configured 
with the parameters in Table 11. Figure 34 shows the BER plots of the system in AWGN 
channel for norma liz ed bandwidth ((3 / ) values from 0.07% to 0.28% in 0.07% incre¬ 

ments. 


#of 

carriers 

#of 

symbols 

Channel 

model 

seed 

Interleaver 

matrix 

a 

range 

p/«, 

<%> 

48 

20000 

1 

33 

fl44 1391 

0-0.8 

0-0.0028 


Table 12. Phase Noise in AWGN Channel Simulation 



Figure 34. BER versus E b /N a Plots for QPSK in AWGN Channel with Norma liz ed 

Bandwidth (3/ R s Values from 0.0% to 0.28%. 
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Figure 34 suggests a linear increase in the required E b /N 0 values to ob¬ 
tain the same error rates as the norma li zed bandwidth (1 / R s was increased. This was an 
expected result because from Equation 3.14, the SNR degradation increases linearly with 
the increasing norma liz ed bandwidth. 

c. Mobile Channel 

This step was performed to investigate the effects of phase noise in a se¬ 
vere mobile channel. The same parameters in Table 12 were used, but the simulation was 
performed under mobile Channel 1. 

Figure 35 shows the BER performance of the system in mobile Channel 1 
for norma li zed bandwidth values ranging from 0 to 0.28% in 0.07% increments. The dot¬ 
ted lines are for the phase noise simulation in the AWGN channel and the solid lines are 
for a phase noise simulation in the mobile Channel 1. From Figure 35, we can observe the 
degradation due to mobile channel by comparing the phase noise simulation results in an 
AWGN channel (dotted lines) with the phase noise simulation results in the mobile chan¬ 
nel (solid lines). The OFDM system under mobile channel conditions requires an E h /N a 
value of approximately 2 dB more than the AWGN case. 
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Figure 35. BER versus E b /N a Plots for QPSK in Mobile Channel 1 with Normalized 

Bandwidth p /R Values from 0.0% to 0.28%. 


C. SUMMARY 

This chapter has described a Matlab-based OFDM communication system that 
can meet the bit rate and performance requirements of the IEEE 802.1 la standard. After 
the verification of the code, several simulation steps were performed under a variety of 
channel conditions. Some of these simulations especially the ones in steps 4 and 5 took 
seven hours, which was a limiting factor in the simulation process. 

The simulation results did not exactly match the theoretical results; however, they 
reflected the expected behavior according to the theory. On the other hand, the results 
were similar to those reported in studies for the same type of simulations. 
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V. CONCLUSION 


The objective of this thesis was to investigate the frequency offset and phase 
noise effects on OFDM based communication systems. This objective was accomplished 
by investigating the effects of frequency offset and phase noise on the performance of the 
system through a comprehensive simulation study. 

A. SUMMARY OF THE WORK DONE 

A comprehensive background of OFDM based communication systems was pre¬ 
sented, including the sub-carrier orthogonality requirement and IFFT/FFT processes. A 
theoretical analysis of frequency errors and their effects on overall system performance 
were examined based on the work reported in the literature. 

A Matlab simulation of an OFDM communication system was developed. The 
simulation studies were performed for the following channels: ideal, AWGN, and mobile. 
The mobile channel includes multipath effects in AWGN. The effects of frequency offset 
and phase noise were studied for all channel models by observing the deviation from the 
AWGN case with increasing frequency offset and phase noise. Simulations studied the 
BER performance plots of OFDM system affected by frequency errors under AWGN and 
mobile channel conditions. 

B. SIGNIFICANT RESULTS AND CONCLUSIONS 

The results were satisfactory in terms of the OFDM system exhibiting the ex¬ 
pected behavior. However, there were some discrepancies with the theoretical results. 
These discrepancies stem from many different factors. To make a one-to-one comparison 
of results from different studies, the simulations must be performed under the same con¬ 
ditions. This was not possible in this thesis because the conditions were either different or 
not known in detail. Therefore, observing the expected trend, rather than obtaining the 
exact numerical results, was enough to realize the objectives of this thesis. 

From the theoretical analyses given in Chapter 3, we have seen that the SNR deg¬ 
radation due to the Doppler shift is generally negligible. The simulations support this 
conclusion. In Figure 26, we can see that mobile Channel 1 and 2 have similar BER 
curves even though the former has 15 Hz and the latter has 5 Hz of Doppler frequency. 
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The theoretical analyses for frequency offset showed that the OFDM system’s 
performance degrades with the square of the relative frequency offset. This result was 
observed in our simulations. From Figure 29, a frequency offset of 0.1% required 1 dB 
more E b /N a value than the AWGN case, and a frequency offset of 0.2% required 2.8 dB 
more E b /N a for the same BER. 

In the phase noise case, the theoretical analyses showed that the OFDM perform¬ 
ance degraded linearly as we increased the normalized bandwidth. The simulation result 
in Figure 34 supports this point. 

Simulation results indicated that mobility causes performance degradation. Fig¬ 
ures 30 and 35 show that the performance of OFDM under mobility conditions is worse 
compared to that of an AWGN channel. 

C. SUGGESTIONS FOR FUTURE STUDIES 

This thesis was conducted under a certain set of conditions. The effects of chang¬ 
ing these conditions on the performance of OFDM systems open a wide range of research 
topics. 

A baseband implementation of the communication system was simulated in this 
thesis. A future study may focus on an actual RF implementation to observe the system’s 
performance in a more realistic scenario. 

Differential encoding/decoding technique was implemented in this thesis. This 
technique may be replaced with coherent detection techniques to study the system’s sen¬ 
sitivity to synchronization issues including the use of pilot tones. This study may help 
evaluate if the added complexity of the coherent techniques pays off in terms of increased 
performance. 

This work focused on the effects of frequency errors and studied only four multi- 
path channel models. A more comprehensive study on the effects of mobile channels may 
be considered in a future work. Such a study will help in the adoption of OFDM for the 
fourth-generation mobile communications. 
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QPSK was used in this thesis. The effects of other signaling schemes on the sys¬ 
tem’s performance may be studied to find an optimum signaling scheme. Also of interest 
would be a study of the effects of changing the OFDM system parameters used in this 
thesis, especially the guard interval and the number of OFDM sub-carriers on the per¬ 
formance of the system. 
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APPENDIX A. MATLAB CODE EXPLANATION 


This appendix explains the function of each m-file and their relations with each 
other. The code is separated into three parts: transmitter, channel and receiver. 

COFDMSIM: This m-file is the main simulation-governing file. The main con¬ 
figuration parameters, such as sub-carrier number, FFT point number, noise variance 
range, multi-path losses, multi-path delays, Doppler frequencies, guard interval length, 
channel models, seed value, interleaver matrix dimensions, time/lfequency differential 
encoding and constellation sizes are chosen here. After the system is configured with 
these parameters, the m-file CHANCDL is called. 

CHANCDL: This m-file performs the simulations depending on the channel 
model chosen in COFDMSIM. This m-file calls CDRCDLFT. 

A. TRANSMITTER 

1. CDRCDLFT 

This is the COFDM encoder with CDL interleaver and is called to compute the 
frequency array of prearranged modulation values and the matrix of differentially en¬ 
coded complex values. It calls MLTPL. 

MLTPL computes the possible interleaver matrix dimensions. CDRCDLFT cre¬ 
ates an error message if the interleaver matrix cannot accommodate all the created sym¬ 
bols. 

MARYMSG generates the required number of bits randomly and then sends them 
to the convolutional encoder. MARYMSG calls CNV_ENCD and BM. 

CNV_ENCD convolutionaly encodes the randomly generated bits in MARYMSG. 
BM then converts the convolutionaly-encoded bits to m-ary symbols. (m value was cho¬ 
sen in COFDM.) 

As the last step, MARYMSG creates the m-ary msg matrix, which contains the 
message symbols. After obtaining the message symbols from MARYMSG, CDRCDLFT 
calls CDLILV, which interleaves the message symbols. During the interleaving opera¬ 
tion, ROTM is called to rotate the input vector. 
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With the creation of the interleaved message matrix, CDRCDLFT calls BM again 
and creates n-ary message matrix that are to be transmitted. (The choice of n was also 
made in COFDM and it determines the constellation. In the simulations, n = 2 for QPSK 
signaling was used.) 

DIFCDRFT performs time or frequency differential encoding, as required. If time 
differential encoding is chosen, a cumulative summation on each column is executed. For 
frequency differential encoding, a cumulative summation on each row is performed. The 
cumulative summation is executed such that the first element is added to the next element 
and replaces the second element in the current summation, then adds the current sum to 
the third element and replaces it also with the current sum. All the sums are modulo-A. 

After deriving the differentially encoded matrix, its elements are encoded as A-ary 
complex modulation values. The unit circle is divided into the number phase partitions 
depending on the signaling constellation and each element of the differentially encoded 
matrix is mapped to one of these phases. A reference row of ones is appended to the mes¬ 
sage array to provide a reference starting point for the receiver, in the case of time differ¬ 
ential encoding, and two reference columns of ones are inserted in the case of frequency 
differential encoding. These reference symbols are stripped off in the receiver during dif¬ 
ferential decoding. 

CMV2FA is the last m-file CDRCDLFT calls. 

CMV2FA rearranges the differentially encoded complex modulation values into a 
frequency array to prepare them for OFDM frequency generation by the IFFT. 

This final step of the encoder block has the convolutionaly encoded, block inter¬ 
leaved, differentially encoded and prearranged frequency array. The frequency array is 
returned to the CHANCDL m-file for the IFFT process. TDA is the m-file responsible for 
the IFFT process. 
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2. TDA 

Performs the IFFT processing and generates the OFDM frequencies. The trans¬ 
mitted OFDM symbols for multi-path effects are prepared by appending a periodic guard 
interval with length Ng. (Ng is chosen in COFDM and, guard interval details were dis¬ 
cussed in Chapter II.) 

Hie time domain samples representing the OFDM symbol appended with the 
guard interval are returned to CHANCDL as the signal to be transmitted through the 
channel. Figure 36 shows the relations of m-files in the transmitter block diagram. 



Figure 36. Hierarchical m-file Block Diagram for the Transmitter 


B. CHANNEL 

One of the three main channel models can be chosen and combined with fre¬ 
quency offset or phase noise. The phase noise can be applied prior to the channel to simu¬ 
late the phase noise case. 

1. Phase Noise Simulations 

The time domain OFDM symbol samples from the TDA m-file are taken and sent 
to the phase noise m-files of PHSNS or phsns2. 

phsns2 generates phase noise as a Wiener process and then applies it to all the 
samples that are being transmitted. Then, CHANCDL sends the time domain samples to 
the desired channel model. Figure 37 illustrates the place of the phase noise block in the 
m-file hierarchy. 
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Figure 37. Phase Noise Creating Block in the m-file Hierarchy 

2. Channel Model Simulations 

Time domain samples are next sent to the channel model of choice. 
a. Channel Model 0 

This is the ideal noise-free channel used for code check simulation pur¬ 
poses. Here, the output of the transmitter block is sent directly to the receiver. This model 
can also be used to check whether frequency offset or phase noise can be implemented. 
These two possibilities are illustrated in Figure 38. 



Figure 38. Channel Model 0 m-file Hierarchy 
b. Channel Model 1 

If channel model 1 is chosen, then CHANCDL directs the time domain 
samples to AWGN. 

AWGN generates a matrix of complex random numbers chosen from a 
normal distribution with mean 0 and variance 1. This matrix has the same dimensions as 
the time domain samples matrix. This matrix is then multiplied by the noise variance pa¬ 
rameter chosen in COFDM and added to the time domain signal matrix. This channel 
model can either be used individually or as a component of a mobile channel. 
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c. Mobile Channel Models 

Four different mobile channels are created by combining an AWGN chan¬ 
nel with multi-path effects. If desired, the multi-path channels can also be used individu¬ 
ally for simulation purposes. The choice of a mobile channel model is determined in the 
COFDM m-ffle. 

Mobile Channel 1 and 2: If one of these mobile channels is chosen in 
COFDM, then CHANCDL sends the time domain signal to CHUHF. CHUHF functions 
as the multi-path component of mobile Channels 1 and 2. It calls D I.IN F to realize a de¬ 
lay line and set up the multiple delayed paths. 

CVDD creates a continuous digital delay element in which the time do¬ 
main samples are filtered using an eight-tap FIR filter whose tap coefficients are a func¬ 
tion of the desired delay. 

After receiving the delayed paths from DLINE, CHUHF next calls OFST, 
which applies a frequency offset as a ratio of the max Doppler shift to the direct path of 
time domain samples. 

RAY_DOP calculates the max Doppler shift frequency as a ratio of fre¬ 
quency spacing. A sequence of complex numbers with zero mean and variance 0.5 is cre¬ 
ated having a Rayleigh envelope with a mean square value of 1. The real and imaginary 
parts are generated independently. 

Finally, CHUHF includes the power losses for the individual multi-paths 
by multiplying each loss amount by the respective delay line output arrays. At this point, 
CHANCDL has the time domain representation of the transmitted signal plus the multi- 
path effects. The time domain samples are then sent to AWGN as the second component 
of the mobile channel. 

Mobile Channel 3: If this model is chosen in COFDM, then CHANCDL 
routes the signal to channel_a, which models different Doppler frequencies, power losses, 
envelope (Ricean) and delays. It calls jakes, which implements FIR filtering by using a 
Hamming window. 
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Mobile Channel 4\ If this model is chosen, CHANCDL sends our signal to 
channel_b, which is the last mobile channel model with very similar characteristics to 
model 3. It also calls jokes for filtering purposes. 

Mobile Channels 3 and 4 also calls AWGN as the second component of 
the mobile channel. 

Figure 39 shows hierarchical architecture of all the mobile channels. For 
each simulation, it is only possible to choose one mobile channel. However, as the white 
blocks suggests, it is possible to combine mobile channels with phase noise or frequency 
offset. Mobile channels can also be used individually. 



Figure 39. Mobile Channels m-file Hierarchy 
Figure 40 summarizes the entire channel model m-file architecture. 


CHANCDL 


TRANSMITTER 


PHSNS2/PHSNS 


Channel Model 


FROFST 


RECEIVER 


Code check 


AWGN 


Mobile Channels 


Figure 40. Summary of Channel Model m-file Architecture 
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3. Frequency Offset Simulations 

If frequency offset simulations are desired, CHANCDL directs the time domain 
samples to frofst, after the channel model of choice. 


frofst applies the frequency offset to each of the transmitted time domain samples 
as a percentage of the frequency spacing. Figure 41 shows the place of frequency offset 
generating block in the m-file hierarchy. 



Figure 41. Frequency Offset Block in the m-file Hierarchy 


Finally, all the m-files in the channel part have been discussed .The last part of the 
code implements the receiver. 

C. RECEIVER 

The first m-file in the receiver part is the ITDA. 

IT DA takes the time domain samples, removes the guard interval and performs 
FFT to create the frequency domain samples. The decoding functions are performed 
within the DECDRCDL. 

DECDRCDL performs the decoding and deinterleaving functions. It calls 
FA2CMA, which reconstmcts the frequency array as a complex modulation array with 
the correct ordering of frequencies.DECDRCDL next calls DFDCDRFT to decode the 
complex modulation values differentially. 

DFDCDRFT takes the complex modulation values from FA2CMA and performs 
inverse mapping of the complex numbers to N- ary symbols based on the value chosen in 
COFDM. Differential decoding is the inverse of encoding performed in the transmitter. 

After receiving A-ary message from DFDCDRFT, DECDRCDL calls MB and 
BM respectively to reformat the received symbols into /n-ary symbols. 
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The next step in the receiver block is to deinterleave the message symbol array. 
DECDRCDL calls CDLDLV, which creates the deinterleaved message matrix by using 
the ROTM as described in the transmitter block. 

These deinterleaved message symbols are next changed to binary and sent to the 
Viterbi decoder. 

VITERBI performs the Viterbi decoding and during the process it calls BIN2DECI 
and DECI2BIN m-files. 


Finally, DECDRCDL takes the received bit sequence from VITERBI, reformats it 
as an m-ary message sequence and creates the received message matrix. Figure 42 shows 
the receiver m-file block diagram. 



Figure 42. Hierarchical m-file Block Diagram for the Receiver 


At this step, CHANCDL successfully governed the transmitter and receiver op¬ 
erations. Further steps are taken to plot the necessary data to be used for analysis pur¬ 
poses. The simulations end after providing the simulations’ performance curves. 


Figure 43 summarizes the entire m-file algorithm. 



Figure 43. General m- file Block Diagram 

66 





































































APPENDIX B. MATLAB CODE 


In this appendix, the m-files used in the simulations are presented. Descriptions of 
the m-files included here can be found in Appendix A. 

%AWGN 

% 

%Title : Additive White Gaussian Noise (Channel Model 1) 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%INPUTS : 

% X : Time domain samples of the transmitted signal 
% s : Seed parameter 

% sigma : AWGN Noise variance parameter for calculating Eb/No 

% 

%OUTPUTS: 

% Y : Time domain samples of the transmitted signal plus AWGN 

% 

function Y = awgn(X, sigma, s) 

% 

%Find dimensions of the input array 
[rr,cc]=size(X); 

% 

randn('seed',s+30); 

%Generate a random real part 
wreal=randn(rr,cc); 

%Generate a random imaginary part 
randn('seed',s+40); 
wimg=i*randn(rr,cc); 

% 

%An array of random complex entries chosen from a normal distribution with 
%mean 0.0 and variance 1.0. Array dimensions are the same as X. 

W=wreal+wimg; 

% 

%Random noise multiplied by the sigma factor and added to the signal. 

Y=X+(sigma. *W); 
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%BIN2DECI 

% 

% 

%Title : Binary To Decimal Conversion 

%Author : Tan Kok Chye, Naval Postgraduate School 

% 

%INPUT : 

% v_x : Binary input 

% 

%OUTPUT: 

% v_y : Decimal output 

% 

function v_y=bin2dcci(v_x) 
v_l=length(v_x); 
v_y=(v_l-l:-l:0); 
v_y=2. A v_y; 
v y=v x*v y'; 
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%BM 

% 

% 

%Title : Binary to M-ary Converter 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 

%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% 

%INPUTS: 

% q : Base 2 exponent for M-ary symbol generation 
% v : Binary data vector 

% 

%OUTPUTS: 

% m : M-ary output vector in decimal notation 

% 

function m=bm(q,v) 

% 

%Find the length of input vector,v,and determine if there is a remainder 

%after dividing by q 

n=length(v); 

r=rcm(n,q); 

% 

%If there is no remainder,don't pad v input vector. Otherwise add the appropriate 
%number of zeros to generate a code word with an exact multiple of q bits. 

% 

if r==0 

v=v; 

else 

v=[v zeros(l,q-r)]; 
end 

% 

%Place least significant bit of the symbol on the left end. 
map=l; 
for j=l:q-l 
map=[map 2 A j]; 
end 
% 

%Remove q bits at a time from v to generate m-ary symbol values. 

n=length(v); 

p=round(n/q); 

A=zeros(q,p); 


A(:)=v; 

m=map*A; 

m_ary_msg=m; 
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%CDLDLV 

% 

% 

%Title : CDL Block Deinterleaver 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% 

%INPUTS: 

% 1 : Intermediate matrix row number 

% k : Intermediate matrix column number 

% dcase: Indicates the deinterleaving method to be 

% used (9 different cases). In our thesis only block interleaving was used. 

% si : Input message string to be deinterleaved 

% SYNC : Frame synchronization bits (Not used in COFDM simulation) 

% 

%OUTPUTS: 

% s : Interleaved output string 

% 

function s=cdldlv(l,k,dcase,si,SYNC) 

si(length(si)+l-length(SYNC):length(si))=zeros(l,length(SYNC)); 

N=length(si); 

if l*k==N 

x=zeros(l,k); 

x(:)=si; 

K=(l:k)-1; 

CR=K.*(K+l)/2; 

L=(l:l)-1; 

RR=L.*(L+l)/2; 

% 

if dcase==l 
for kk=l:k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

elseif dcase==2 
for kk=l:k 

[z,x(: ,kk)] =rotm(x(: ,kk) ,CR(kk)); 
end 

elseif dcase==3 
for kk=l:l 

x(kk,: )=rotm(x(kk,: ),RR(kk)); 
end 

elseif dcase==4 
for kk=l:l 

[z,x(kk,:)] =rotm(x(kk,:) ,RR(kk)); 
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end 

elseif dcase==5 
for kk=l:k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

for 11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

elseif dcase==6 
for kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

for 11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

elseif dcase==7 
for kk=l:k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

for 11=1:1 

[z,x(ll,:)J=rotm(x(ll,:),RR(ll)); 

end 

elseif dcase==8 
for kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

for 11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

end 

x=x'; 

s=x(0; 

s=s'; 

end 
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%CDLILV 

% 

% 

%Title : CDL Block Interleaver 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% 

%INPUTS: 

% 1 : Intermediate matrix row number 

% k : Intermediate matrix column number 

% dcase: Indicates the deinterleaving method to be 

% used (9 different cases). In our thesis only block interleaving was used. 

% si : Input message string to be deinterleaved 

% SYNC : Frame synchronization bits (Not used in COFDM simulation) 

% 

%OUTPUTS: 

% si: Interleaved output string 

% 

%Subroutines Used: rotm.m 

% 

function si = cdlilv(l,k,dcase,s,SYNC) 

N=length(s); 

if l*k==N 

x=zeros(l,k); 

x=x'; 

x(:)=s; 

x=x'; 

Intermediate_mx=x; 

K=(l:k)-1; 

CR=K.*(K+l)/2; 

L=(l:l)-1; 

RR=L.*(L+l)/2; 

% 

if dcase==l 
for kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

elseif dcase==2 
for kk= 1 :k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

elseif dcase==3 
for kk=l:l 

[z,x(kk,:)] =rotm(x(kk,:) ,RR(kk)); 
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end 

elseif dcase==4 
for kk=l:l 

x(kk,: )=rotm(x(kk,: ),RR(kk)); 
end 

elseif dcase==5 
for 11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

for kk=l:k 

[z,x(:,ldc)]=rotm(x(:,kk),CR(kk)); 

end 

elseif dcase==6 
for 11=1:1 

[z,x(ll,:)]=rotm(x(ll,:),RR(ll)); 

end 

for kk=l:k 

x(:,ldc)=rotm(x(:,kk),CR(kk)); 

end 

elseif dcase==7 
for 11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

for kk=l:k 

[z,x(:,kk)]=rotm(x(:,kk),CR(kk)); 

end 

elseif dcase==8 
for 11=1:1 

x(ll,:)=rotm(x(ll,:),RR(ll)); 

end 

for kk=l:k 

x(:,kk)=rotm(x(:,kk),CR(kk)); 

end 

end 

x; 

si=x(:); 

si=si'; 

end 

si(length(si)-length(SYNC)+l:length(si))=SYNC; 
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%CDRCDLFT 

% 

% 

%Title : COFDM Encoder with CDL Interleaver 
%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% 

%INPUTS: 

% pic : Parameter indicating the figure number 
% s : Seed parameter 

% freqno : Number of OFDM frequencies (sub-carriers) used in each message array 
% rintlv : Parameter for intermediate matrix row number 
% cintlv : Parameter for intermediate matrix column number 
% N : Number of EFT frequency sample points,must be larger than freqno 
% mary : Initial M-ary symbol format (OFDM symbol bit length) 

% nary : Final N-ary symbol format (PSK symbol bit length) 

% fort : Selects either frequency (fort=l) or time (fort=0) differential encoding 
% 

%OUTPUTS: 

% Fa : Frequency array of prearranged modulation values 
% MD : Matrix of differentially encoded complex values 
%(unity magnitude and one of N-ary possible phases (N-PSK)) 

% B : Matrix of 8-ary symbols 

% nsymno: Number of N-ary generated symbols 

% 

%Subroutines Used: marymsgm,cdlilv.m,mb.m,bmm,difcdrft.m,cmv2fa.m 

% 

func- 

tion[Fa,MD,B_ce,B_random,nsymno]=cdrcdlft(pic,dcase,s,freqno,rintlv ,cintlv,N,mary,na 
ry,fort); 

% 

% Determine whether the number of OFDM frequencies are even, indicated 
% by the "freqno" parameter. If odd, go to error message. Odd frequencies are not al¬ 
lowed 

% since the formation of the frequency array is symmetrical and even. 

% 

if rem(freqno,2)~=0 

disp(ERROR: The number of matrix columns, representing OFDM frequencies, must be 
an even number!') 
elseif rem(freqno,2)==0 

% 

% Check if interleaver matrix dimensions are greater than freqno. 

% If not, then display error message and stop. 

% 
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if (rintlv*cintlv)<(freqno) 
disp(") 

disp('ERROR: The row and column interleave parameters are not compatible with # of 
OFDM frequencies!') 
disp(") 
else 

% Calculate the row symbol number 
symno=rintiv*cintiv/freqno; 

% 

% Display error message if symno and freqno not compatible with rintlv and cintlv and 
stop. 

% If not compatible,the interleaver function does not work correctly. 

% 

if rem(symno, 1) ~=0 
disp(") 

disp('ERROR: The row and column interleave parameters are not compatible with # 
of OFDM frequencies!') 

disp(' For the entered rintlv, cintlv, and freqno parameters, the calculated symno is:') 
disp(symno) 

multiesall=mltpl(rindv,cintlv); 
multies=multiesall( 1 ,(2:length(multiesall) -1)); 
disp(' Possible choices for freqno based upon rintlv and cintlv are:') 
disp(") 
disp(multies) 
elseif rem(symno,l)==0 
if freqno >= N; 
disp(") 

disp('ERROR: The number of frequency points, N, needs to be increased !') 

disp('N must be larger than:') 

disp(") 

disp(freqno) 

disp(") 

elseif freqno < N; 

% 

% Generate a random message matrix of m-ary symbols,based upon parameter,mary. 
Nmbr_of_symbols =symno*freqno; 

[B_ce,B_random]=marymsg(mary,symno,freqno); 
Rndm_m_ary_msg=B_random; 

% 

% Perform a block interleaving function on the matrix, B, with rintlv rows 
% and cintlv columns. 

SYNC=[]; 

[Br Bc]=size(B_ce); 

Bt=B_ce'; 

Bvect=Bt(:)'; 
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si=cdlilv(rintlv,cintlv,dcase,Bvect,SYNC); 

Bi=reshape(si,Bc,Br)'; 

Lrtrlvd_array=Bi; 
ml =bm(nary ,mb(mary ,Bi)); 
lengthml=length(ml); 
nsymno=lengthml; 
remml=rem(lengthml ,freqno); 
if remml=0; 
ml=ml; 
else 

zero=zeros(freqno-remml); 

ml=[ml zero(l,:)]; 

end 

Iength2ml=length(ml); 

m=(reshape(ml,freqno,length2ml/freqno))'; 

N_ary_msg=m; 

% 

% Generate a differentially encoded matrix of complex 
% values with unity magnitude and one of (2 A n) equal phases. 
MDD=difcdrft(nary ,m,fort); 

[MDm MDn]=size(MDD); 

MD=MDD; 

Cmplx_mod_array=MDD; 

% 

% Form the frequency array of modulation values that include guard interval. 
Fa=cmv2fa(N,MD); 

Freq_array=Fa; 

end 

end 

end 

end 
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%CHANCDL 

% 

% 

%Title : Simulations for AWGN & Multipath Fading Channel 
%Author : Dave Roderick, Naval Postgraduate School 
%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%Subroutines Used: cdrcdlft.m,tda.m,awgn.m,chuhf.m,itda.m,decdrcdl.m,check.m 

% 

function 

chancdl(chnmdl,wait,pic,dcase,s,freqno,rintlv,cintlv,N,mary,nary ,n,k,blklgth,Ng,sigs,loss 

,dly,dop,freqspace,fort) 

sigvect=sigs; 

klgth=length(k); 

chklp=l; 

symno=rintlv*cintlv/freqno; 

%%%%FOR PHASE NOISE%%%%% 

%for perctB=0: .0008: .0032 

%%%%FOR FREQUENCY OFFSET%%%% 

%for donci-0:0.0007: .0028 

errvect=[]; 

bervect=[]; 

freqerrmx=[]; 

errsperpr=[]; 

Es_No=[]; 

Eb_No=[]; 
sermx=[]; 
bermx=[]; 
rowerrmx=[]; 

for lp=l:length(sigvect); 

[xmt,modvals,B_ce,B_random,nsymno]=cdrcdlft(pic,dcase,s,freqno,rintlv,cintlv,N,mary, 
nary,fort); 

xmtifft=tda(Ng,xmt); 

%%%%%%%%%%%%FOR PHASE NOISE%%%% 
%sandnaphns=phsns2(xmtifft,freqspace,perctB); 

%xmtifft=sandnaphns; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

xmtpts=l:length(xmtifft); 
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if chnmdl=0 
sandn=xmtifft; 

elseif chmndl=l 

disp(['Sigma=',num2str(sigvect(lp))]); 

sandn=awgn(xmtifft,sigvect(lp),s); 

elseif chnmdl=21 

sandrnltpth=chuM(s+l,xmtifft,loss,dly,dop,N,freqspace); 

disp(['Sigma=',num2str(sigvect(lp))]); 

sandn=awgn(sandmltpth,sigvect(lp),s); 

elseif chnmdl=22 

sandmltpth=chuhf(s+l ,xmtifft,loss,dly,dop.N,freqspacc); 

disp(['Sigma=',num2str(sigvect(lp))]); 

sandn=awgn(sandmltpth,sigvect(lp),s); 

elseif chnmdl=23 

disp(['Sigma=',num2str(sigvect(lp))]); 
sandmltpth=channel_a(xmtifft); 
sandn=awgn(sandmltpth,sigvect(lp),s); 

elseif chnmdl=24 

disp(['Sigma=',num2str(sigvect(lp))]); 
sandm 1 tpth=chan ncl_b( x mt ifft); 
sandn=awgn(sandmltpth,sigvect(lp),s); 
end 

%%%%%%%%%FOR FREQUENCY OFFSET%%%%%% 

%sandnofst=frofst(sandn,doner); 

%sandn=sandnofst; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

sandnfft=itda(Ng. sandn); 

K=(length(modvals( 1,: )))/2; 

[rcvd,rcvd_hit,random_msg,random_hit.M.MM]=dccdrcdl(pic,dcasc,K,sandnfft,nsymno, 
freqno,rintlv,cindv,rnary,nary,fort ,B_random); 

%%Transmitted_msg=B_random; 

Transmitted_msg=random_msg; 

Received_msg=rcvd; 

[er¬ 
rors,bit_error,freqerrs,errmx,rowerrs]=check(pic,random_msg,random_bit,rcvd,rcvd_bit, 

n,k(chklp),blklgth); 


78 



errvect=[errvect,errors]; 
bervect=[bervect,bit_error]; 
freqerrmx=[freqerrmx;freqerrs]; 
rowerrmx=[rowerrmx;rowerrs]; 

cmtEs_No=l/(2*N*(sigvect(lp) A 2)); 
cmtEb_N o=cmtEs_N o; 

% %based on M=4 i.e. for coded QPSK, Eb=Es. 

Es_No=[Es_No,cmtEs_No]; 

Eb_No=[Eb_No,cmtEb_No]; 

Es_Nodb= 10*log 10(Es_N o); 

Eb_Nodb=10*log 10(Eb_No); 
end 

sei-en'vcct/(symno*frcqno); 
ber=bervect/(2*symno*freqno); 
sermx=[sermx;ser]; 
bermx=[bermx;ber]; 
errsum=sum(errvect); 
errsperpr=[errsperpr,errsum]; 
errmax=max(rowerrmx'); 

% 

% PLOTS 

% 

figure(pic+l) 
plofrmodvals,'*') 
hold on; 
plot(0,0,'+') 
hold off; 

%title(['Transmitted Signal Constellation',int2str(2 A nary), '-ary (QPSK)']) 

xlabel('Re'); 

ylabel('Im'); 

axis('square');orient tall;grid 
pause(wait); 

%%%%Plot of FFT POINTS%%%% 

%figure(pic+2) 

%plot([0:N - l],abs(xmt),'*') 

%title(['Frequency Array Plot (number of FFT frequency points are ',int2str(N),')']) 
%xlabel(['Guard interval length is ',int2str(N-freqno)]) 

%axis('square');orient tall;grid 
%pause(wait); 

% 

figure(pic+2) 
surf(abs(modvals)); 
shading interp;grid;orient tall 
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%title(['Magnitude Plot of Transmitted Signal']) 
xlabel('OFDM Subcanier Number') 
ylabel('Symbol Row Number') 
zlabel('Signal Magnitude') 
pause(wait); 

% 

figure(pic+3) 
plot(M,'*');hold on; 
plot(0,0,'+');hold off; 

%title(['Received Signal Constellation', int2str(2 A nary),'-ary (QPSK), (Before Differential 

Decoding)']) 

xlabel('Re'); 

ylabel('Im'); 

orient tall;axis('square');grid 
pause(wait); 

% 

figure(pic44) 
plot(MM,'+') 
hold on; 
plot(0,0,V) 
hold off; 

%title(['Received Signal Constellation',int2str(2 A nary),'-ary (QPSK), (After Differential 

Decoding)']) 

xlabel('Re'); 

ylabel('Im'); 

orient tall;axis('square');grid 
pause(wait); 

% 

figure(pic+5) 
surf(abs(M)); 
shading interp 
grid;orient tall 

%title(['Magnitude Variation of Received Signal (Sigma=',num2str(sigvect(lp)),')']) 

xlabel('OFDM Subcanier Number') 

ylabel('Symbol Row Number') 

zlabel('Signal Magnitude') 

pause(wait); 

% 

if errsum~=0 

%%%%PERFORMANCE PLOTS%%%% 

%theoretical on top of simulation curves 
%q=length(Eb_No) 

%for i=l:q 

%P=0.5*(erfc(Eb_No(l,i))*(l-0.25*erfc(Eb_No(l,i)))); 
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%Pb(l,i)=P; 

%end 

figure(pic+6) 

%semilogy(Eb_Nodb,Pb,'m');grid;hold on; 
semilogy(Eb_Nodb,ber);hold on;grid 
if fort==l 

%tide('Performance curve: BER vs. Eb/No for QPSK Signaling (Freq. Diff.Enc.)') 
elseif fort==0 

%tide('Performance curve: BER vs. Eb/No for QPSK Signaling (Time. DiffEnc.)') 
end 

xlabel('Eb/No(dB)'); 

ylabel('BER'); 

%legend('theorical curvcfsimulation curve');orient tall 

% 

figure(pic+7) 

%semilogy(Es_Nodb,2*Pb,'m');hold on 
semilogy(Es_Nodb,ser);grid;hold on 
if fort==l 

%title('Performance curve: SER vs. Es/No (Freq. Diff.Enc.)') 
elseif fort==0 

%title('Performance curve: SER vs. Es/No (Time. Diff.Enc.)') 
end 

xlabel('Es/No(dB)'); 

ylabel('SER'); 

%legend('theoretical curve','simulation curve') ;orient tall 
end 
end 
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%channel_a 

% 

% 

%Title : Mobile Channel 2 

%Author : Prof.Roberto Cristi, Naval Postgraduate School 
%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%Subroutines Used: jakes.m 

% 

function yr=channcl_a(xt) 

[asd asn]=size(xt); 

xt=xt.'; 

xt=xt(:).'; 

% 

% yr=channel(xt) 

% with xt=transmitted vector 
% yr=received vector 

% length(yr)=length(xt) 

% 

% if the program gives you an error, just increase the size of the 
% transmitted vector "xt" 

Fs=l .25* 10 A 6; % sampling frequency (Hz) 

Fd=200; % doppler frequency (Hz) 

%Parameters for C 
Td=[0.0, 0.25,0.5, 

PdB=[0.0,-1.0,-9. 

K=[0,0,0,0,0,0]; 

Td=Td*(10 A (-6)); 
nd=round(Td*Fs); 

P=10. A (PdB/10); 

Np=length(xt); 

yr=zcros(sizc(xt)); 

for k= 1 :length(Td) 
g=jakes(Fd, Fs, Np); 

s=sqrt(P(k)/(K(k)+1)); % random path 

m=sqrt(P(k)*K(k)/(K(k)+l)); % direct path 
total=s*g+m*exp(j *2*pi*Fd/Fs*(0:length(xt)-1)); 
yr=yr+[ zcrost 1 ,nd(k)), xt( 1 :length(xt)- nd(k))]. *total; 
end 

yr=reshape(yr,asn,asd).'; 


% time delays (in microsec) 
% Powers (in dB) 

% Ricean Factor 

% time delays (in sec) 

% time delay in samples 

% Powers (Finear) 


Channel A 
,1.0,1.9,2.2]; 

0,-10.0,-15.0,-20.0]; 
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%channel_b 

% 

% 

%Title : Severe Mobile channel 

%Author : Prof.Roberto Cristi, Naval Postgraduate School 
%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%Subroutines Used: jakes.m 

% 

function yr=channcl_b(xt) 

[asd asn]=size(xt); 

xt=xt.'; 

xt=xt(:).'; 

% 

% yr=channel(xt) 

% with xt=transmitted vector 
% yr=received vector 

% length(yr)=length(xt) 

% 

% if the program gives you an error, just increase the size of the 
% transmitted vector "xt" 

Fs=l .25* 10 A 6; % sampling frequency (Hz) 

Fd=200; % doppler frequency (Hz) 

% Parameters for Channel B 

Td=[0.0, 0.25,9.0,13.0, 17.0, 20.0]; % time delays (in microsec) 

PdB=[-2.5,0.0,-12.8,-10, -25.2, -16]; % Powers (in dB) 
K=[0.5,0.5,0,0,0,0]; % Ricean Factors 

Td=Td*(10 A (-6)); % time delays (in sec) 

nd=round(Td*Fs); % time delay in samples 

P=10. A (PdB/10); % Powers (Finear) 

Np=length(xt); 

yr=zcros(sizc(xt)); 

for k= 1 :length(Td) 
g=jakes(Fd, Fs, Np); 

s=sqrt(P(k)/(K(k)+1)); % random path 

m=sqrt(P(k)*K(k)/(K(k)-i-l)); % direct path 
total=s*g+m*exp(j *2*pi*Fd/Fs*(0:length(xt)-1)); 
yr=yr+[ zeros! 1 ,nd(k)), xt( 1 :length(xt)- nd(k))]. *total; 
end 

yr=reshape(yr,asn,asd).'; 
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%CHECK 

% 

% 

%Title : Source and Sink Message Checker 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 

%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% - 

func- 

tion[error_no,bit_error_total,freqerrs,errmx,rowerrs]=check(pic,x,xbit,y,ybit,n,k,blklgth) 

if blklgth>n 

disp(") 

disp(ERROR! The block length(blklgth), must be equal to or less than the code word 
length(n)') 

disp('Please enter a smaller value for blklgth, or change n') 
disp(") 

elseif blklgth<=n 
if n<k 
disp(") 

disp('Error! The code word length(n) must be equal to or larger than the information 
length(k)') 

disp('Please enter a larger value for n, or change k to a smaller number.') 

disp(") 

elseif n>=k 

First_matrix=x; 

Second_matrix=y; 

[rx cx]=size(x); 

% 

%Compare inputs x and y and generate error matrix, "errors" 
errors=(x~=y); 

First=xbit; 

Second=ybit; 

[rxl cxl]=size(xbit); 

% 

%Compare inputs xbit and ybit and generate BIT error matrix, "bit_errors" 
bit_errors=(xbit~=ybit); 

% 

%Find the error distribution vs. OFDM frequencies 
freqerrs=sum(errors); 

% 

%Find the error location in "errors" where element in x and y differ. 
Error_locations=(find(errors))'; 

Errorjiumber=sum(sum(errors)); 

Correct_symbl_num=(size(y, 1 )*size(y,2))- Error_number; 

% 
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%Find the bit error location in "errors" where element in xl and yl differ. 
%bit_Error_locations=(find(bit_errors))'; 

bit_error_total=sum(sum(bit_errors)); 

%Correct_bit_num=(size(y 1,1 )*size(y 1,2))-bit_Error_number; 

%Reed-Solomon 8-bit symbol correction for (n-k)/2 symbols 
symcon-floor) (n- k)/2); 
if blklgth<=(n-k) 

disp('Error!! !The block length is too short for the given n and k values') 
disp(") 

elseif blklgth>(n-k) 
errtrans=errors'; 

% 

%Reshape the error matrix as a vector of errors 
errvect=errtrans(:; 
blkrem=rem(length(errvect),blklgth); 
if blkrem~=0; 

zeropad=zeros(blklgth- blkrem); 
errvectpad=[errvect zeropad(l,:)]; 
elseif blkrem==0; 
errvectpad=errvect; 
end 
% 

blknos=length(errvectpad)/blklgth; 

errcorct=[]; 

errblksum=[]; 

% 

for lp=l:blknos; 

errblk=errvectpad(((blklgth*(lp -1))+1): (blklgth*lp)); 
errblklgth=length(errblk); 

if sum(eiTblk)<=symcorr; 
noen-zeros) crrblklgth); 
errblk=noerr(l,:); 

elseif sum(en’blk)>symcorr; 
errblk=errblk; 
end 

errcorct=[errcorct errblk]; 
eiTblksum=[errblksum sum(errblk)]; 
end 

newerrvect=errcorct( 1 :length(errvect)); 
errtot=sum(newerrvect); 

RSerrs=(reshape(newerrvect,size(errors,2),size(errors,l)))'; 

% 

%Find the error distribution vs. OFDM Frequencies 
freqerrs=sum(RSerrs); 
enindex=(find(RSerrs))'; 
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RSerrtot=sum(errblksum); 

RSerrdif=Error_number-RSerrtot; 
errperblk=[( 1 :blknos) ;errblksum]; 

% 

%Check if x and y are the same. If not, display error message 
if x=y; 

disp('There are no errors!!!') 

error_no=0; 

errmx=errors; 

rowerrs=sum(errors'); 

else 

disp('WARNING!:Errors were detected!') 

disp(") 

if n=k 

disp('WARNING!: Since n=k,there is no R-S error correcting possible') 

disp(") 

end 

disp(['For the given input parameters:n=',int2str(n),'and k=',int2str(k),',the Reed- 
Solomon code is capable']) 

disp(['of correcting ',int2str(symcorr),'errors.']) 
disp(") 

% 

%Check if xbit and ybit are the same. If not, display error message 
if xbit==ybit; 

disp('There are no bit errors!!!') 
bit_error_no=0; 
bit_errmx=bit_errors; 
bit_rowerrs=sum(bit_errors'); 
else 

disp('WARNING!:Errors were detected!') 
disp(") 
if n=k 

disp('WARNING!: Since n=k,there is no R-S error correcting possible') 
disp(") 
end 
end 

% 

%RS code was able to correct all errors 
if errtot==0 

Pre_RS_error_matrix=errors; 

disp('EXCELLENT: Tlie Reed-Solomon code corrected all detected errors!') 

disp(['Originally the error total was:',int2str(Error_number)]) 

disp(") 

error_no=0; 

errmx=zeros(rx,cx); 
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ro werrs=sum(errmx'); 


% 

%RS code was able to correct some errors but not all of them 
elseif errtot<Error_number 
Pre_RS_error_matrix=errors; 

Post_RS_error_matrix=RSerrs; 
errmx=RSerrs; 
rowerrs=sum(errmx'); 

disp('The Reed-Solemon code corrected some detected errors, but not all.') 

disp(['Originally the error total was : ',int2str(Error_number)]) 

disp(") 

disp(['After R-S decoding , the error number was reduced to:',int2str(RSerrtot)]) 
disp(") 

error_no=RSerrtot; 

disp(['The total number of correct symbols are:',int2str((size(y,l)*size(y,2))-RSerrtot)]) 
disp(") 

disp('The error number distribution per block number is :') 
disp(errperblk) 

% 

%RS code did not correct any errors 
elseif errtot==Error_number 
Error_matrix=errors; 
errmx=errors; 
rowerrs=sum(errors'); 

disp('The Reed-Solomon code did not correct any errors.') 
disp('Perhaps a more powerful R-S code is required.') 
disp(") 

disp(['The total number of error occurrences is:',int2str(Error_number)]) 
disp(") 

error_no=errtot; 

disp('The error number distribution per block number is :') 
disp(errperblk) 
end 
end 
end 
end 
end 

disp('.__'); 
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%CHUHF 

% - 

%Title : UHF Channel Model (multipath Channel Model 1) 
%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

function y=chuhf(s,x,loss,dly,dop,N,freqspace) 

c=10. A (-loss./20); 

deltat= l/(N*freqspace); 

d=(dly. * .000001) ./deltat; 

e=dop ./freqspace; 

[L,Nt]=size(x); 

D=length(d); 

x=x.'; 

x=x(:).'; 

% 

%D path with delays from d. 
xd=dline(x,d); 

[rr,cc]=size(xd); 

x=xd(l,:); 

% 

% Offsets direct path by .7 of max doppler freq. 
xo=ofst( .7 *e( 1) ,N,x); 

% 

% First path with no fading, 
for 1=1 :D 

a=ray_dop(s ,cc,N,e(l)); 

xd(l,:)=a.*xd(l,:); 

end 

%Sums the fading paths 
y=c*xd; 

% 

%Adds in the First path without fading 

y=y+xo; 

y=y(l:F*Nt); 

y=reshape(y,Nt,F).'; 
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%CMV2FA 

% 

% 

%Title : Complex Frequency Array Generator 
%Author : Dave Roderick, Naval Postgraduate School 
%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%Changes complex modulation values to frequency array 

% 

function X=cmv2fa(N,M) 

[m n]=size(M); 
if rem(n,2)==0; 

M=M; 

else 

% 

M=[zeros(m, 1) M]; 
end 

[m n]=size(M); 

K=round(n/2); 

% 

%Generate a matrix of zeros with m row and N columns. 
X=zeros(m,N); 

% 

X(:,1:K)=M(:,K+1:2*K); 

X(:,N-K+1:N)=M(:,1:K); 
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%CNV_ENCD 

% 

% 

%Title : Convolutional Encoding 

%Reference : Contemporary Communication System using Matlab 
%John G. Proakis & Masoud Salehi. 

% 

function ce_output=cnv_encd(ce_g,ce_kO,ce_input) 

% cnv_encd(ce_g,ce_kO,ce_input) 

% determines the output sequence of a binary convolutional encoder 
% ce_g is the generator matrix of the convolutional code 
% with ce_nO rows and ce_l*ce_kO columns. Its rows are ce_gl,ce_g2,....ce_gn. 
% ce_kO is the number of bits entering the encoder at each clock cycle. 

% check to see if extra zero padding is necessary 
if rem(length(ce_input),ce_kO)>0 

ce_input=[ce_input,zeros(size( 1 :ce_kO-rem(length(ce_input),ce_kO)))]; 
end 

ce_n=length(ce_input)/ce_kO; 

%check the size of matrix ce_g 
if rem(size(ce_g,2),ce_k0)>0 
error('Error, ce_g is not of the right size.') 
end 

% determine ce_l and ce_nO 
ce_l=size(ce_g,2)/ce_k0; 

%disp(['The value of ce_l is:',int2str(ce_l)]); 
ce_nO=size(ce_g, 1); 

%disp('') 

%disp(['The value of ce_nO is:',int2str(ce_n0)]); 

%add extra zeros 

ce_u=[zeros(size( 1 :(ce_l-1 )*ce_kO)),ce_input,zeros(size( 1 :(ce_l-1 )*ce_kO))]; 

%generate ce_uu, a matrix whose column are the contents of 

%conv. encoder at various cycles. 

ce_ul=ce_u(ce_l*ce_kO:-1:1); 

for ce_i= 1 :ce_n+ce_l- 2 

ce_ul=[ce_ul,ce_u((ce_i+ce_l)*ce_kO:-l:ce_i*ce_kO+l)]; 

end 

ce_uu=reshape(ce_u 1 ,ce_l*ce_kO,ce_n+ce_l-1); 

%determine the ce_output 

ce_output=reshape(rem(ce_g*ce_uu,2),l,ce_n0*(ce_l+ce_n-l)); 
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clc,close all 
%COFDMSIM 

% 

% 

%Title : Simulation Of COFDM 

%Author : Dave Roderick, Naval Postgraduate School 

%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% 

disp('___'); 

disp('This batch m-file mns COFDM simulations using different channel models.') 
%fort=input('To run the frequency version, enter l(one), To mn the time version, enter 
O(zero), or to mn both enter 2(two):'); 
fort=l; %frequency version 

%freqno=input('Enter the # of OFDM frequencies (note : must be even):'); 
freqno=48; 

%%%%N=input('Enter the number of EFT points (Note : This number must be larger 
than # of OFDM frequencies):'); 

N=64; 

chnmdl=input('Choose the channel model; 0-(Noise Free), l-(AWGN), 21-(Mobile 
Channel-1), 22-(MC-2), 23-(MC-3), 24-(MC-4)? :'); 
if chnmdl=0 

disp('Code Check simulation.'); 
sigs=0;loss=0;dop=0;dly=0; 
elseif chnmdl=l 

disp('AWGN Channel simulation.'); 

sigs=input('Enter the noise variance(sigma) range or single value. (Ex lin- 
space(0,0.02,20)or .003):'); 

loss=0;dop=0;dly=0; 
elseif chnmdl=21 

disp('Mobile Channel-1.'); 

sigs=input('Enter the noise variance(sigma) range or single value. (Ex lin- 
space(0,0.02,20)or .003):'); 

loss=[0,2.17,4.34,6.51,8.69,10.86,13.03,15.20,17.37,19.54,21.71,23.89,26.06,28.23,30.4, 
32.57,34.74,36.92]; 

dop=[15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15,15]; 

dly=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0. 
85]; 

elseif chnmdl=22 

disp('Mobile Channel-2.'); 

sigs=input('Enter the noise variance(sigma) range or single value. (Ex lin- 
space(0,0.02,20)or .003):'); 
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loss=[0,2.17,4.34,6.51,8.69,10.86,13.03,15.20,17.37,19.54,21.71,23.89,26.06,28.23,30.4, 
32.57,34.74,36.92]; 

dop=[5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5]; 

dly=[0,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65,0.70,0.75,0.80,0. 
85]; 

elseif chnmdl==23 
disp('Mobile Channel-3.'); 

sigs=input('Enter the noise variance(sigma) range or single value. (Ex lin- 
space(0,0.02,20)or .003):'); 

loss=0;dop=0;dly=0; 
elseif chnmdl=24 

disp('Mobile Channel-4.'); 

sigs=input('Enter the noise variance!sigma) range (Ex: linspace(0,0.02,20)or .003):'); 
loss=0;dop=0;dly=0; 
end 

%%%% Choosing Interleaver %%%% 

%allcase=input('Simulate all interleaver cases (yes) or specific ones(no)? (l=yes,0=no):'); 

allcase=0; 

if allcase=l 

disp('All cases,(0-8),will be tested.'); 
cases=[0:8]; 
elseif allcase==0 

%%%%cases=input('Enter specific case numbers from (0 to 8)(Ex [0 4 5 8]):'); 

cases=0; % Block interleaving chosen 
end 


if fort~=2 

if length(cases)~=l 

casey_n=input('Do you want to find optimal interleaver case(s) ? (l=yes, 0=no):'); 
else 

casey_n=0; 

end 

end 

%%%%% 

totsym=input('Enter the total minimum number of symbols to simulate (Ex 10000):'); 
rowno=ceil(totsym/freqno); 
if totsym~=(rowno*freqno) 

disp(['Note:Tlie actual total number of symbol to be simulated will be 
:',int2str(rowno*freqno)]); 
end 

%pry_n=input('For the interleaver, do you want to calculate all possible intermediate ma¬ 
trix dimension pairs?(l=yes,0=no):'); 
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pry_n= 1; 
pairl=l; 

pair2=rowno*freqno; 
if pry_n== 1 

Intrlvr_pairs=intlvprs(rowno,ffeqno); 
intlvrprs=Intrl vr_pairs; 

%%%% Desired interleaver pairs can be entered %%%% 

%disp(") 

disp('For these input parameters, all possible interleaver dimension pairs are:') 
disp(Intrlvr_pairs) 
end 

pairs=input([Desired interleaver pair? (Ex [row # col #] = [20 50] (Note: enter- 
ing[',int2str(pairl),' ',int2str(pair2),'],or [',int2str(pair2),' ',int2str(pairl),'], offers no inter¬ 
leaving functionality):']); 

%rintlv=intlvrprs(8,l);cintlv=intlvrprs(8,2) 

rindv=pairs(l); 

cindv=pairs(2); 

%mary=input('Enter the number of M-ary bits, q (i.e. for 256-ary, q=8):'); 
mary=l; 

%nary=input('Enter the number of N-ary bits,q(i.e. for 16-ary, q=4):'); 
nary=2; % QPSK is chosen 
%freqspace=round( 16600000/freqno); 

freqspace=312500; % According to the standard 0.3125 Mhz of Frequency spacing. 
%Ng=input('Enter the guard interval length (Number of sample points):'); 

Ng=16; 

%ecc=input('Do you want to include Reed Solomon error correction coding ? (l=yes, 
0=no):'); 

ecc=0;% Reed Solomon Error correction coding is not chosen 
if ecc==l 

%code=input('Enter n,k and error correction block length (Ex [240 200 240]):'); 
code=[240 200 240]; 
n=code(l); 
k=code(2); 
blklgth=code(3); 
elseif ecc==0 
n=freqno; 
k=freqno; 
blklgth=freqno; 
end 

svals=input('Enter specific seed values, or 0 for a random seed (ex [103 22, 60] or [0]):'); 

wait=3;pic=0; 

svect=[]; 

for mn=l:length(svals); 
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errvect=[]; 
errcase=[]; 
errtot=[]; 
if min(svals)=0 

rand('seed',sum(100*clock)); 
s=round(abs(rand( 1 )*pi* 10*(pic+1 )*run)); 
elseif min(svals)~=0 
s=svals(run); 
end 

svect=[svect,s]; 
for 1=1 iength(cases); 

dispC.__') 

disp(['Run#:',int2str(run)]); 
disp(['Seed=',int2str(s)]); 
disp( ['Interleaver case=',int2str(cases(l))]); 
if fort<=l 


chancdl(chnmdl,wait,pic,cases(l),s,freqno,rintlv,cintlv,N,mary,nary ,n,k,blklgth,Ng,sigs,lo 
ss,dly,dop,freqspace,fort); 
elseif fort==2 

disp('Frequency differential encoding/decoding simulation...') 

disp('.__') 


chancdl(chnmdl,wait,pic,cases(l),s,freqno,rintlv,cintlv,N,mary,nary ,n,k,blklgth,Ng,sigs,lo 
ss,dly,dop,freqspace, 1); 

disp(************************************************************) 

disp('Time differential encoding/dccoding simulation....') 

dispC_.') 


chancdl(chnmdl,wait,pic+12,cases(l),s,freqno,rintlv,cintlv,N,mary,nary ,n,k,blklgth,Ng,sig 
s,loss,dly,dop,freqspaee,0); 
end 
end 
end 

) P( 

disp(");disp('Simulation finished!') 

Seed=svect; 
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%CVDD 

% 

% 

%Title : Continuous Variable digital delay element. 

%Reference : C.W. Farrow," A Continuously Variable Digital Element", IEEE 
%Intemational Symposium on Circuits & Systems,pp.2641-2645,1988. 

% 

% . 

function [y]=cvdd(x,alpha) 
if ((nargin~=2)l(nargout~=l)) 
error('ERROR:usage: y=y=cvdd(x,alpha);'); 
return; 
end 

if (size(x)~=size(alpha)) 

error('ERROR:x and alpha must be the same size'); 

return; 

end 

if (abs(alpha)>0.5) 

errorf'ERROR:alpha must be within -0.5 and 0.5'); 

return; 

end 

% 

%—----- 

% In itia liz ation 

% 

% 

% In itia li ze FIR filter coefficients are in [1] (0,0.328 pass band) 

C0=[-0.013824 0.054062 -0.157959 0.616394 0.616394 -0.157959 0.054062 -0.013824]; 
Cl=[0.003143 -0.019287 0.1008 -1.226364 1.226364 -0.1008 0.019287 -0.003143]; 
C2=[0.055298 -0.216248 0.631836 -0.465576 -0.465576 0.631836 -0.216248 0.055298]; 
C3=[-0.012573 0.077148 -0.403198 0.905457 -0.905457 0.403198 -0.077148 0.012573]; 

% 

% -- - .... 

% 4 parallel FIR and add together based on [1] 

% 

y0=filter(C0, [ 1 ] ,x); 
yl=filter(Cl,[l],x); 
y2=filter(C2,[l],x); 
y3=filter(C3, [ 1 ] ,x); 

% 

y=alpha.*y3; 

y=alpha.*(y+y2); 

y=alpha.*(y+yl); 

y=y+y0; 
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%DECDRCDL 

% 

% 

%Title : COFDM Decoder With Deinterleaveing 
%Author : Dave Roderick, Naval Postgraduate School 
%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%-— 

func- 

tion[outmsg,viterbi_output_bit,random_msg,random_bit,M,MM]=decdrcdl(pic,dcase,K,F 

a,nsymno.freqno,rdintlv,cdintlv.mary.nai'y,foit.B_random) 

% ’ 

M=fa2cma(K,Fa); 

Cmplx_mod_vals=M; 

% 

naryp=nary; 

[s,MM]=dfdcdrft(naryp,nary ,M,fort); 

[L,cc]=size(s); 
strans=s'; 
svect=strans(:).'; 
corrs=s vect( 1: nsy mno); 

% 

nsymno; 

Br=bm(mary,mb(nary ,corrs)); 
lengthBi-length(Br); 
rmndn=rem(length(Br),freqno); 
if rmndi^=0; 

Bi-Br; 

elseif rmndr~=0; 

Bi-Bif 1 :(lengthBr- rmndr)); 
end 

rcvd=(reshape(Br,freqno,length(Br)/ffeqno))'; 

Rcvd_Intlv_Ary=rcvd; 

% 

[Br Bc]=size(rcvd); 

SYNC=[]; 

sr=rcvd'; 

si=sr(:)'; 

sd=cdldlv(rdintlv,cdintlv,dcase,si,SYNC); 
received=reshape(sd,Bc,Br)'; 
viter_G=[l 0 1 1 0 1 1 ;1 1 1 1 0 0 1]; 
vitcr_k=l; 

binary_value=mb(mary ,sd); 

[viterbi_output,survivor_sta,cumul_metrix]=viterbi(viter_G,viter_k,binary_value); 
mary_dec=bm(mary ,viterbi_output); 
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viterbi_output_bit=viterbi_output; 

%outmsg=reshape(sd,Bc,Br)'; 

% 

random_bit=B_random; 
random_msg=bm(mary ,random_bit); 
[Brow Bcol]=size(random_msg); 

% 

outmsg=reshape(mary_dec,Bcol,Brow)'; 

Sink_Msg=outmsg; 
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%DECI2BIN 

% 

% 

%Title : Decimal to Binary Converter 

%Author : Tan Kok Chye, Naval Postgraduate School 

% 

function y=deci2bin(x,l) 

y=zeros(l,l); 

vi=l; 

while x>=0 & vi<=l 

y(vi)=rem(x,2); 

x=(x-y(vi))/2; 

vi=vi+l; 

end 

y=y(i:- i:i); 
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%DFDCDRFT 

% 

% 

%Title : Complex Number Demodulator & Frequency/Time Differential Decoder 
%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% 

function [s,M]=dfdcdrft(qp,q,MD,fort) 
if fort==0 %Time Differential decoding 

% 

MD=MD'; 

[m n]=size(MD); 

% 

% Perform a looping routine to find the phase differences between adjacent values in the 
% array,MD,and put these calculated values into array,M. 
for 1=1 :m 
for j=l:n-l 

M(l,j)=MD(l,j+l)*conj(MD(l,j)); 

end 

end 

% 

%Transpose the array back to its original form 
M=M'; 

% 

% Calculate the number of M-ary symbols based upon the exponent qp,then use this 
number 

% to find the number of equally spaced phases in a unit circle. 

N=2 A qp; 

dph=2*pi/N; 

% 

% Divide the phase arguments of elements in M, by the equal phases generated by dph. 
phn=angle(M) ./dph; 

% 

% Calculate the phase sector number by finding the remainders. 

s=rem(round(phn)+N,N); 

elseif fort==l % Frequency Differential decoding 

% 

% Transpose the modulation array, and find the dimensions 
[m,n]=size(MD); 

MD=MD(:,2:n); 

[m n]=size(MD); 

% 

% Perform a looping routine to find the phase differences between 
% adjacent values in the array, MD, and put these calculated values into array,M. 
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for 1=1 :m 
for j=l:n-l 

M(l,j )=MD(l,j+1 )*conj (MD(l,j)); 
end 
end 
N=2 A qp; 
dph=2*pi/N; 

% 

% Calculate the phase sector number by finding the remainders. 

phn=angle(M) ./dph; 

s=rem(round(phn)+N,N); 

end 
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%DIFCDRFT 

% 

% 

%Title : Complex Number Modulator & Frequency/Time Differential Encoder 
%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

function MD=difcdrft(q,m,fort) 
if fort==0 %Time differential encoding 
% M-ary alphabet size 
N=2 A q; 

% Determine the number of equal phases based upon the m-ary symbol length 
dph=2*pi/N; 

% Find the size of the input symbol matrix 
[rr n]=size(m); 

% 

% Perform the time differential encoding of phase values by cumulative summation, 

% down one column at a time across the entire matrix. This function generates a matrix, 
for k=l:n 

md=cumsum( m( :,k)); 

% Generate the complex numbers with corresponding phase values. 

MD(: ,k)=exp(i*dph.*md); 
end 

% 

% Inject the reference row of ones (zero phase) at top of output matrix for 
% differential encoding synchronization 
MD=[ones(l,n); MD]; 

elseif fort==l % Frequency Differential encoding 

% 

% M-ary alphabet size 

N=2 A q; 

dph=2*pi/N; 

% Find the size of the input symbol matrix 
[rr n]=size(m); 

% 

md=cumsum(m'); 

md=md'; 

% 

% Generate the complex numbers with corresponding phase values. 

MD=exp(i*dph. *md); 

% 

% Inject the reference row of ones (zero phase) at top of output matrix for 
% differential encoding synchronization. 

MD=[ones(rr,2) MD]; 
end 
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%DIFFCHKR 

% 

% 

%Title : Differential Encoder/Decoder Checker 
%Author : Dave Roderick, Naval Postgraduate School 
%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% - 

function diffchkr(s,symno,freqno,mary,nary) 

fort=input(Eor the frequency version, enter 1 (one); for the time version, enter 0 (zero):'); 

% 

B=marymsg(mary,s,symno,freqno); 

Rndm_m_ary_msg=B; 

% 

ml=bm(nary,mb(mary,B)); 

lengthml=length(ml); 

m=(reshape(ml,lengthml/symno,symno))'; 

N_ary_msg=m; 

% 

if fort==l 

disp(");disp('Frequency Differential Encoding/Decoding version') 

% 

%Freq. Diff. Enc. 

% 

MDD=difcdrf(mary,m); 
elseif fort~=l 

disp(");disp('Time Differential Encoding/Decoding version') 

% 

MDD=difcdrt(mary ,m); 
end 

% 

maryq=mary; 
if fort==l 

% 

[s M]=difdcdrf(maryq,mary,MDD); 
elseif fort ~=1 

% 

[s M]=difdcdrt(maryq,mary,MDD); 
end 

% 

%Check results for correctness. 

[error_no,freqerrs,errmx,rowerrs]=check(0,m,s,freqno,freqno,freqno); 
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%DLINE 

% 

% 

%Title : UHF Channel Delay L ine Generator 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 

% 

function xd=dline(x,d) 
x=x.'; 

dmax=max(d); 
dmin=min(d); 
nmm=floor(drnin); 
nmax=ceil(dmax); 
x=[x;zeros(nmax+3,1)]; 

N=length(x); 

Nd=length(d); 

% 

for n=l:Nd; 
di=d(n); 

D=floor(di); 

dcld=di-D; 

xd(:,n)=cvdd(x,deld- .5); 
xd(:,n)=[zeros(D,l);xd(l :N-D,n)]; 
end 

xd=xd.'; 

[rr,cc]=size(xd); 
xd=xd(: ,4+nmin:cc); 
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%FA2CMA 

% 

% 

%Title : Frequency Array To Complex Modulation Array Converter 
%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%—-----...—. 

function Mm=fa2cma(K,X) 

[m n]=size(X); 

Mm(:,l:K)=X(:,n-K+l:n); 

Mm(: ,K+1:2*K)=X(:, 1 :K); 

Cmplx_mod_vals=Mm; 
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%FROFST 

% 

%Title : Frequency Offset 

%Author : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% 

function [xofstd]=frofst(x,y) 

[sati sutu]=size(x); 
x=x.'; 
x=x(:).'; 
wo=length(x); 
wu=l:wo; 
yd=wu*y 

xofst=x. *exp(i*(2.5 *pi). *yd); 
xofst=xofst( 1: sutu* sati); 
xofst=reshape(xofst,sutu,sati).'; 
xofstd=xofst; 
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%INTLVCHK 

% 

% 

%Title : Interleaver/Deinterleaver Verifier 

%Author : Dave Roderick, Naval Postgraduate School 

%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% - 

function intlvchk(s,synmo,freqno,rintlv,cintlv,mary,dcase) 
multiples=mltpl(symno,freqno); 

Intrlvr_nbr_mltpls=multiples; 

% 

if (symno*freqno)~=(rintlv*cintlv) 

disp('ERROR: The interleaver matrix dimensions are not compatible with the message 
array size.') 

disp(' Possible matrix dimensions are:');disp(") 
dispfmultiples) 

disp('Note: The selected matrix dimensions can not accomodate the message array.') 
disp(' In this case the number of rows times the number of columns is:');disp(") 
disp(symno*freqno) 
elseif(symno*fieqno)/(rintlv*cintlv)= 1 

% 

B=maiymsg(mary,s,symno,freqno); 

Random_msg=B 

% 

SYNC=[]; 

[Br Bc]=size(B); 

Bt=B'; 

Bvect=Bt(:)'; 

si=cdhlv(rintlv,cintlv,dcase,Bvect,SYNC); 

Bi=reshape(si,Bc,Br)'; 

Interleaved_array=Bi 

% 

[Br Bc]=size(Bi); 

SYNC=[]; 

sr=Bi'; 

si=sr(:)'; 

sd=cdldlv(rintlv,cintlv,dcase,si,SYNC); 

Bd=reshape(sd,Bc,Br)'; 

Deinterleaved_array=Bd 

% 

[error_no,fi'eqerrs,errmx,rowerrs]=check(0,B,Bd,freqno,freqno,freqno); 

end 
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%INTLVPRS 

% 

% 

%Title : Intermediate Matrix Interleaver Dimension Pairs 
%Author : Dave Roderick, Naval Postgraduate School 
%Modified By : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%-— .-.—.—---- 

function pairs=intlvprs(n,m) 
prod=n*m; 
multvect=[l]; 
for i=2:prod; 

remdr=rem(prod,i); 
if remdi-=0 

multvect=[multvect ij; 
else 

multvect=multvect; 

end 

mult=multvect; 

end 

lngth=length(mult); 
nhi-mulUlngth); 
result=[l nbr]; 
for i=2:lngth; 

cmtpr=[mult(i) nbr/mult(i)]; 
result=[result;cmtpr]; 
end 

pairs=result; 
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%ITDA 

% 

% 

%Title : Frequency Domain Samples Without Guard Interval 
%Author : Prof. Paul H. Moose, Naval Postgraduate School 

% 

% 

function Y=itda(Ng,y) 

% 

[L Nt]=size(y); 

% Remove the guard interval for channel compensation. 

% 

y=y(:,Ng+l:Nt); 

% Take the FFT of array, y 

% 

Y=fft(y.').'; 
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%JAKES 

% 

% 

%Title : Multi-path 

%Author : Prof. Roberto Cristi, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%—- 

function g=jakes(Fd,Fs,N); 

% g=jakes(Fd,Fs,N) 

% generate a random vector of length N with jakes spectrum 
% Fd=doppler frequency 
% Fs=sampling frequency (Fs>10*Fd required) 

% N=vector length 

% g=generated random vector 

% The vector g has unit power, ie. g'*g*Fs/N=l 


F=round(log2(Fs/( 10*Fd))); 

F0=Fs/(2 A F); % samphng frequency for generating the random sequence 

% 

N0=2*ceil(N/(2 A F)); 

Nf=512; % number of frequency components in spectrum definition 

Nh=256; % length of FER filter 

% FIR Filter impulse response 

w=0.54-0.46*cos(2*pi*(0:Nh- 1)/Nh); % hamming window 

fd=Fd/F0; % digital doppler frequency 

f=0:(l/Nf):l-(l/Nf); % vector of digital frequencies (l/2=Nyquist Freq.) 

kd=floor(fd*Nf); % index for doppler frequency fd 

H0( 1 :kd)=sqrt( 1 - (f( 1 :kd)/fd). A 2); 

I=find(H0==0); H0(I)=0.0001*ones(size(I)); 

H=zeros(l,Nf); 

H( 1 :length(H0))=( 1 ./HO); 

H=H.*exp(-j*2*pi*f*Nh/2); 
hO=real(ifft(H)); 
h=hO( 1 :length(w)).*w; 
end 

% Generate time varying taps at low sampling freq Fs_ch 
seed=14; 

randn('state',seed); 
a=randn( 1 ,N0+Nh); 
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seed=seed+6; 
b=randn( 1 ,NO+Nh); 

%x=randn( 1 ,NO+Nh)+j *randn( 1 ,NO+Nh); 
x=a+j*b; 

g=filter(h,l,x); 

g=g(Nh+l:Nh+NO); % steady state response 


% Upsample to Fs=(2 A L)*F0 in L stages 

for m=l:L 

F0=2*F0; 

g=reshape([g;zeros(size(g))], l,2*length(g)); 
omegad=2*pi*Fd/F0; 

Domega=(pi/2)- omegad; 
omegac=((pi/2)+omegad)/2; 

M=ceil(((8*pi/Domega)-1 )/2); 

Nfilt=2*M+1; % filter order using hamming window 

nt=0:Nfilt-1; 

w=0.54-0.46*cos(2*pi*nt/Nfilt); % hamming window 

hm( 1 :M)=sin(omegac*(nt( 1 :M)- M)) ./(pi*(nt( 1 :M)- M)); 
hm(M+1 )=omegac/pi; 

hm(M+2:Nfilt)=sin(omegac*(nt(M+2:Nfilt)-M))./(pi*(nt(M+2:Nfilt)-M)); 

hm=hm(l :Nfilt).*w; 

g=filter(hm,l,g); 

g=g(Nfilt:length(g)); % steady state 
end 

% Normalize 
g=g(l:N); 

K=g*g7N; 

g=g/sqrt(K); 
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%MARYMSG 

% 

% 

%Title : M-ary Message Test Pattern Generator 
%Author : Tan Kok Chye, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% - 

function [vmary_ce4'andom_bit]=marymsg(q,n,m) 

% 

[random_bit]=msg(n*m*q*.5-6); 

conv_g=[l 0 1 101 1;1 1 1 1001]; 
conv_k0=l; 

conv_output=cnv_encd(conv_g,conv_kO,random_bit); 
decml=bm(q,conv_output); 
deciml=decml( 1,1: (n*m)); 
vmary_ce=(reshape(deciml,m,n))'; 
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%MB 

% 

% 

%Title : M-ary To Binary Converter 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 

% 

function [b]=mb(q,m) 

% 

row=size(m,l); 

col=size(m,2); 

% 

m=reshape(m', 1 ,(row*col)); 

% 

b0=rem(m,2); 

m=(m-b0)./2; 

B=bO; 

% 

for j=l:q-l 

bj=rem(m,2); 

m=(m-bj)./2; 

B=[B;bj]; 

end 

% 

b=B(:)'; 

binary=b; 
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%METRIC 

%-.—.-----. 

% 

%Title : Viterbi Hard Decision Decoding metric 
%Author : Tan Kok Chye, Naval Postgraduate School 

% 

function distance=metric(v_x,v_y) 

if v_x==v_y 

distance=0; 

else 

distance=l; 

end 
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%MLTPL 

% 

% 

%Title : Common Multiples 

%Author : Tan Kok Chye, Naval Postgraduate School 

% 

function [mult]=mltpl(n,m) 

max=n*m; 

multvect=[l]; 

% 

for i=2:max; 
remdr=rem(max,i); 
if remdi-=0 
multvect=[multvect ij; 
else 

multvect=multvect; 

end 

mult=multvect; 

end 
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%MSG 

% 

% 

%Title : Message Test Pattern Generator 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 

%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% - 

function u=msg(k) 

% 

%rand('uniform'); 

% 

%temp=rand('seed'); 

% 

%rand('seed',s); 

% 

u=randint(l,k); 
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%NXT_STAT 

% 

% 

%Title : Next State 

%Author : Tan Kok Chye, Naval Postgraduate School 

% 

function [next_state,memory_contents]=nxt_stat(current_state,input,v_L,v_k) 
binary_state=deci2bin(current_state, v_k*(v_L-1)); 
binary_input=deci2bin(input,v_k); 

next_state_binary=[binary_input,binary_state(l:(v_L-2)*v_k)]; 

next_state=bin2deci(next_state_binary); 

memory_contents=[binary_input,binary_state]; 
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%OFST 

% 

% 

%Title : Channel Offset 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% - 

function xo=ofst(e,N,x) 

[m Nt]=size(x); 

xo=x.'; 

x=x.'; 

x=x(:); 

x=x.'; 

Nt=length(x); 

1=1 :Nt; 

%Creating the offset frequency 

% 

ex=x.*exp(i*(2*pi/N)*e.*l); 

xo(:)=x; 

xo=xo.'; 
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%phsns2 

% 

% 

%Title : Phase Noise 

%Author : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%This m.file creates phase noise as a Wiener Process for all time samples 

% 

function phsnsig=phsns2(xmtifft,freqspace,perctB) 
perctB; 

[str stn]=size(xmtifft); 

T =str*stn; 

B=freqspace*perctB; 
randn('state', 100) ; 
dW = zeros(l,stn*str); 

W = zeros(l,stn*str); 

% 

dW(l) = normrnd(0,sqrt(2*B)); 

W(l) = dW(l); 
for j = 2:(stn*str) 
dW(j) = normrnd(0,sqrt(2*B*j)); 

W(j) = (W(j-1) + dW(j)); 
end 

xmtifft=xmtifft.'; 

xmtifft=xmtifft(:).'; 

xmtifftpn=xmtifft.*exp(i*(2*pi*4* 10 A (- 6)). *W); 

xmtifftpn=reshape(xmtifftpn,stn,str).'; 

phsnsig=xmtifftpn; 

% phase noise plot 
figure( 14) 

plot([0:T],[0,W],'r-') 

xlabel('f) 

ylabel('PN(t)') 

title('Phase Noise for all Time Domain Samples') 
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%phsns 

% 

%Title : Phase Noise 

%Author : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

%This m.file creates phase noise as a Wiener Process for each OFDM symbol 

% 

function phsnsig=phsns(xmtifft,freqspace,perctB) 

[str stn]=size(xmtifft); 

T =stn; 

B=freqspace*perctB; 

for i=l:str 

randn('state', 100) ; 
dW = zeros(l,stn); 

W = zeros(l,stn); 

dW(l) =normrnd(0,sqrt(2*B)); 

W(l) = dW(l); 
for j = 2:stn 

dW(j) = normrnd(0,sqrt(2*B*j)); 

W(j) = (W(j-1) + dW(j)); 
end 

xmtifftpn(i,:)=xmtifft(i,:).*exp(i*(2*pi*4*10 A (-6)).*W); 

end 

phsnsig=xmtilftpn; 

% phase noise plot 
figure( 14) 

plot([0:T],[0,W],'r-') 

xlabel('t') 

ylabel('PN(t)') 

title('Phase Noise Created for Each OFDM Symbol') 
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%RAY_DOP 

% 

% 

%Title : Rayleigh Doppler 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% -------- 

function c=ray_dop(s,M,N,es) 

m=0:M-l; 

randn('seed',s+10); 

prl=randn(l,20); 

randn('seed',s+20); 

pim=i*randn( 1,20); 

p=prl+pim; 

P= P /(40 a .5); 
rand('seed',s+30); 
e=rand(l,20); 
e=es*cos(2*pi*(e- .5)); 

E=exp(i*2*pi*e' *m/N); 
c=p*E; 
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%ROTM 

% 

% 

%Title : Rotate Vector 

%Author : Prof. Paul H. Moose, Naval Postgraduate School 

% 

function [vp,vn]=rotm(v,m) 

L=length(v); 

m=rem(m,L); 

ii=(l:L)-l; 

isp=rem(ii- m+L,L)+1; 
isn=rem(ii+m+L,L)+l; 
vp=v(isp); 
vn=v(isn); 
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%TDA 

% 

% 

%Title : Time Domain Samples With Guard Interval Precursor 
%Author : Prof. Paul H. Moose, Naval Postgraduate School 
%Modified by : Ahmet Yasin ERDOGAN, Naval Postgraduate School 

% 

% - 

function x=tda(Ng,X) 

[m N]=size(X); 

% 

% Perform inverse EFT on frequency values in array,X 
x=ifft(X.'); 

% Add precursor of Ng samples to the beginning of the time domain array for channel 
% compensation. 

% 

x=x.'; 
if Ng=0 
x=x; 
else 

x=[x(:,N-Ng+l:N) x]; 
end 
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%VITERBI 

% 

% 

%Title : Viterbi Decoder 

%Reference : Contemporary Communication System using MatLab 
%John G. Proakis & Masoud Salehi. 

% 

func- 

tion[decoder_output,survivor_state,cumulated_metric]=viterbi(v_G,v_k,channel_output) 

% 

v_n=size(v_G,l); 

% check the sizes 
if rem(size(v_G,2),v_k)~=0 
error('Size of v_G and v_k do not agree') 
end 
% 

if rem(size(channel_output,2),v_n)~=0 
error('channel output not of the right size') 
end 

v_L=size(v_G,2)/v_k; 
number_of_states=2 A ((v_L-1) * v_k); 

%generate state transition matrix, output matrix, and input matrix 
for vj=0:number_of_states-1 
for v_l=0:2 A v_k-1 

[next_state,memory_contents] =nxt_stat(v _j, v_l,v_L,v_k); 
input(v j+1 ,next_state+1 )=v_l; 
branch_output=rem(memory_contents*v_G',2); 
nextstate(v _j+1, v_l+1 )=next_state; 
output(v j+1 ,v_l+1 )=bin2deci(branch_output); 
end 
end 

state_metric=zeros(number_of_states,2); 

depth_of_trelhs=length(channel_output)/v_n; 

channel_output_matrix=reshape(channel_output,v_n,depth_of_trenis); 

survivor_state=zeros(number_of_states,depth_of_trelhs+l); 

%start decoding of non-tail channel outputs 
for v_i= 1 :depth_of_trellis - v_L+1 
flag=zeros( 1 ,number_of_states); 
if v_i<=v_L 

step=2 A ((v_L- v_i)*v_k); 

else 

step=l; 

end 

for vj=0:step:number_of_states-1 
for v_l=0:2 A v_k-1 
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branch_metric=0; 

binary_output=deci2bin( output( v J+1, v_l+1), v_n); 
for v_ll= 1 :v_n 

branch_metric=branch_metric+metric(channel_output_matrix(v_ll,v_i),binary_output(v_ 

m 

end 

if((state_metric(nextstate(v _j+1 ,v_l+1)+1,2)>state_metric(v J+1,1 )+branch_metric)lflag( 
nextstate(v J+1, v_l+1)+1 )==0) 

state_metric(nextstate(v _j+1 ,v_l+1)+1,2)=state_metric(v J+1,1 )+branch_metric; 
survivor_state(nextstate( v J+1 ,v_l+1)+1 ,v_i+1 )=vJ; 
flag(nextstate(v _j+1 ,v_l+1)+1 )= 1; 
end 
end 
end 

state_metric=state_metric(:,2:-1:1); 
end 

%start decoding of the tail channel_outputs 
for v_i=depth_of_trellis - v_L+2 :depth_of_trellis 
flag=zeros( 1 ,number_of_states); 

last_stop=numbcr_of_statcs/(2 A ((v_i-depth_of_trcllis+v_L-2)*v_k)); 
for vj=0:last_stop-1 
braneh_metric=0; 

binary_output=deci2bin(output(v j+1,1), v_n); 
for v_ll=l:v_n 

branch_metric=branch_metric+metric(channel_output_matrix(v_ll,v_i),binary_output(v_ 

ii)); 

end 

if((state_metric(nextstate(v J+1,1)+1,2)>state_metric(v j+1,1 )+branch_metric) lflag(next 
state(vj+l,l)+l)==0) 

state_metric(nextstate(v _j+1,1)+1,2)=state_metric(v _j+1,1 )+branch_metric; 
survivor_state(nextstate(v J+1,1)+1 ,v_i+1 )=v J; 
flag(nextstate(v J+1,1)+1 )= 1; 
end 
end 

state_metric=state_metric(:,2:-1:1); 
end 

%generate the decoder output from the optimal path 
state_sequence=zeros( 1 ,depth_of_trellis+1); 

state_sequence( 1 ,depth_of_trellis)=survivor_state( 1 ,depth_of_trelhs+1); 
for v_i=l:depth_of_trellis 
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state_sequence( 1 ,depth_ofJrellis - 

v_i+l)=survivor_state((state_sequence( 1 ,depth_of_trellis+2- v_i)+1) ,dcpthoftrcl 1 i s - 

v_i+2); 

end 

decoder_output_matrix=zeros(v_k,depth_of_trellis - v_L+1); 
for v_i=l:depth_of_trellis-v_L+l 

dec_output_deci=input(state_sequence( 1 ,v_i)+1, state_sequence( 1, v_i+1)+1); 
dec_output_bin=deci2bin(dec_output_deci,v_k); 
decoder_output_matrix(:,v_i)=dec_output_bin(v_k:-1:1)'; 
end 

decoder_output=reshape(decoder_output_matrix, 1 ,v_k*(depth_of_trellis - v_L+1)); 
cumulated_metric=state_metric( 1,1); 
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